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Introduction 


Radio  frequency  discharges  currently  play  a  major  role  in  the 
microelectronics  industry  for  the  fabrication  of  Integrated  circuits.  These 
discharges  are  typically  used  for  stripping  of  photoresist,  anisotropic 
etching  of  semiconductor,  oxide,  and  metal  surfaces,  and  the  deposition  of 
organic  and  inorganic  thin  films.  In  all  of  these  applications,  it  is  beneficial 
to  have  a  good  understanding  of  the  physical  and  chemical  processes 
involved  in  these  discharges  and  at  the  surfaces  exposed  to  the  discharge. 
While  there  is  a  long  history  of  extensive  research  on  amd  application  of  dc 
discharges,  only  recently  have  rf  discharges  become  of  significant  interest. 
This  report  describes  theoretical  and  experimental  studies  undertaken  in  an 
effort  to  develop  a  better  understanding  of  some  of  the  underlying  gas 
discharge  ph3rsics  occurring  in  these  discharges. 

Section  II  of  this  report  presents  an  equivalent  circuit  model  for  a  planar 
rf  plasma  reactor.  The  ph3rsical  properties  of  the  plasma  sheaths  adjacent  to 
the  electrodes  are  incorporated  in  the  model.  The  sheath  capacitances  and 
the  conduction  currents  through  the  sheaths  predicted  by  the  model  are 
time  varying  and  have  a  highly  nonlinear  dependence  on  the  potentials 
across  the  plasma  sheaths.  Numerical  calculations  based  on  this  model  show 
that  that  the  waveforms  of  the  voltage  differences  across  the  sheaths  are 
highly  nonsinusoidal  and  agree  with  reported  measurements.  The  results 
given  this  section  were  reported  at  the  37th  Annual  Gaseous  Electronics 
Conference  in  Boulder,  Colorado,  October,  1984  and  have  been  published  in 
Physical  Review  A. 

Section  III  outlines  a  model  for  the  energy  distribution  of  ions 
bombarding  electrode  surfaces  in  planar  rf  plasma  reactors  at  low  pressures, 
where  collisions  in  the  sheath  may  be  ignored.  Given  the  waveform  for  the 
voltage  difference  across  the  sheath  adjacent  to  an  electrode,  the  model 
predicts  the  ion  bombardment  energy  distribution  in  the  limit  where  the 
trzinsit  time  of  the  ions  across  the  sheath  is  much  shorter  than  the  period  of 
the  applied  'f  voltage  and  in  the  limit  where  the  ion  transit  time  is  much 
longer  than  the  rf  period.  Numerical  calculations  utilizing  the  sheath  voltage 
waveforms  from  the  equivalent  circuit  model  presented  in  Sec.  II  show  that 
the  ion  bombardment  energy  distribution  is  blmodal  and  skewed  toward 
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lower  ion  bombardment  energies.  The  results  given  in  this  section  were 
also  presented  at  the  37th  Annual  Gaseous  Electronics  Conference  in 
Boulder,  Colorado.  October,  1984,  and  are  being  submitted  for  publication  to 
the  Journal  of  Vacuum  Science  and  Technology. 

In  an  effort  to  obtain  a  reliable  experimental  data  base  for  comparison 
with  the  results  from  the  ion  bombardment  energy  distribution  model,  an 
experimental  system  was  designed  and  constructed  to  measure  the  ion 
bombardment  energy  distribution  at  an  electrode  in  a  planar  rf  reactor.  The 
design  of  and  resulting  data  from  this  experimental  system  is  reported  in 
Sec.  rv.  This  system  was  designed  so  as  to  fulfill  the  required  assumptions 
of  ion  bombardment  energy  distribution  model.  During  the  design  and 
operation  of  this  system,  several  unanticipated  problems  arose  which  were 
of  Interest  in  themselves  and  are  therefore  discussed.  Section  IV  also 
includes  results  of  measurements  made  with  the  final  functioning  design 
and  discusses  several  conclusions  based  on  these  results. 

In  developing  the  models  presented  in  Secs.  II  and  III,  the  plasma  sheath 
model  used  incorporated  the  Bohm  crlterlum,  which  is  only  applicable  to 
plasmas  consisting  of  a  positive  ion  species,  electrons,  and  neutrals. 

However,  plasmas  containing  negative  ion  species,  in  addition  to  these  other 
species,  are  the  plasmas  of  primary  interest  to  the  microelectronics 
Industry  due  to  the  use  of  electronegative  gases  in  most  plasma  assisted 
processing  steps.  In  an  effort  to  investigate  the  effects  of  negative  ions  on 
the  plasma  sheath  characteristics,  a  sheath  model  was  developed  for  a 
plasma  containing  a  negative  ion  species,  a  positive  ion  species,  electrons, 
and  neutrals.  This  model  is  discussed  in  Sec.  V.  Under  appropriate 
assumptions,  a  criterion  is  derived  which  must  be  satisfied  in  order  for  the 
model  to  be  physically  realistic.  This  criterion  is  then  used  to  calculate  the 
conductivity  and  capacitance  of  a  plasma  sheath  in  the  presence  of  negative 
ions.  The  implications  of  these  results  for  dc  and  rf  discharges  are 
investigated,  including  an  explanation  of  the  confinement  of  negative  ions 
inside  a  plasma  volume.  This  material  was  presented  at  the  39th  Annual 
Gaseous  Electronics  Conference  in  Monterey,  California,  October,  1986,  and 
is  being  submitted  for  publication  to  Physical  Review  A. 

Finally,  a  model  Is  presented  in  Sec.  VI  for  the  radial  distribution  of  ion 
and  electron  densities  and  electric  field  strength  in  the  positive  column  of  a 
dc  discharge  for  a  plasma  consisting  of  a  positive  ion  species,  electrons,  and 
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neutrals..  The  set  of  equations  involved  consists  of  the  particle  and 
momentum  conservation  equations  for  the  ions  and  electrons  and  Poisson's 
equation.  Utilizing  this  single  set  of  equations  and  appropriate  assumptions, 
this  model  has  been  solved,  through  suitable  numerical  techniques.  The 
results  from  this  work  show  the  development  of  both  the  amblpolar  region 
In  the  bulk  of  the  positive  column  and  the  sheath  region  near  the  discharge 
wall.  The  importance  of  the  application  of  this  model  to  discharges  at 
various  pressures  and  tube  radii  is  discussed.  This  work  was  presented  at 
the  39th  Annual  Gaseous  Electronics  Conference  in  Monterey,  California. 
1986,  and  is  being  submitted  for  publication  to  Physics  of  Fluids. 
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Application  of  the  physics  of  plasma  sheaths  to  the  modeling  of 
rf  plasma  reactors 

A.  Metze D.  W.  Ernie,  and  H.  J.  Oskam 

Department  of  Electrical  Engineering,  University  of  Minnesota.  Minneapolis,  Minnesota  55455 
(Received  16  June  1986;  accepted  for  publication  22  July  1986) 

An  equivalent  circuit  model  is  presented  for  a  planar  rf  plasma  reactor.  The  physical  properties 
of  the  plasma  sheath  adjacent  to  the  electrodes  are  incor|>orated  in  the  model.  The  sheath 
capacitances  and  the  conduction  currents  through  the  sheaths  are  time  varying  and  have  a 
highly  nonlinear  dependence  on  the  potentials  across  the  plasma  sheaths.  The  model  shows 
that  the  waveforms  of  the  voltage  differences  across  the  sheaths  are  highly  nonsinusoidal  and 
agree  with  reported  measurements. 


I.  INTRODUCTION 

The  application  of  rf-excited  gaseous  discharges  in  thin- 
film  fabrication  technology  has  become  standard  duri''g  the 
past  decade.  The  simplest  geometry  most  commonly  used  is 
that  of  two  planar  electrodes  between  which  the  rf  voltage  is 
applied. '  A  schematic  representation  of  such  a  planar  rf  plas¬ 
ma  reactor  is  shown  in  Fig.  1.  The  plasma  is  separated  from 
each  electrode  by  a  plasma  sheath.  The  origin  of  these  plas¬ 
ma  sheaths  has  been  discussed  by  several  authors.^-* 

Positive  ions  produced  in  the  plasma  volume  are  accel¬ 
erated  across  the  plasma  sheaths  and  arrive  at  the  electrodes 
with  an  energy  distribution  which  is  determined  by  the  mag¬ 
nitude  and  the  waveform  of  the  time-dependent  potential 
difference  across  the  sheaths,  the  gas  pressure,  etc.  This  ion 
bombardment  energy  distribution  determines  the  degree  of 
anisotropy  in  ihin-film  etching,  amount  of  ion  impact  in¬ 
duced  damage  to  surfaces,  etc.  Therefore,  it  is  important  to 
obtain  a  better  understanding  of  the  dependence  of  the  ion 
bombardment  energy  distribution  on  parameters  such  as 
reactor  geometry,  frequency  of  the  applied  voltage,  gas  pres¬ 
sure,  etc.  Suzuki  et  al*  and  Tsui’  have  reported  models  relat¬ 
ing  to  the  energy  distribution  of  positive  ions  arriving  at  the 
electrodes  of  a  rf  plasma  reactor.  These  authors  assumed 
that  the  electric  potential  difference  across  the  plasma 
sheath  is  simply  the  superposition  of  a  dc  bias  and  a  sinusoi¬ 
dal  voltage. 

In  order  to  obtain  a  more  accurate  model  for  the  ion 
bombardment  energy  distribution,  it  is  necessary  to  develop 
a  model  which  predicts  the  waveforms  of  the  voltages  across 
the  plasma  sheaths  in  a  planar  rf  plasma  reactor.  In  an  effort 
to  determine  these  waveforms,  several  authors  have  used  a 
combination  of  circuit  elements  for  the  representation  of  the 
plasma  sheaths  in  equivalent  circuit  models  of  a  planar  plas¬ 
ma  reactor.  Koenig  and  Maissel*’  used  a  static  capacitor  in 
parallel  with  a  diode.  The  capacitor  represented  the  sheath 
capacitance,  while  the  diode  approximated  the  effect  of  elec¬ 
trons  reaching  the  electrodes.  Horwitz^  added  a  parallel  re¬ 
sistor  to  the  previous  model,  while  other  authors  have  used  a 
simple  parallel  combination  consisting  of  a  resistor  and  ca¬ 
pacitor  of  fixed  values.*  ’  The  latter  model  is  sometimes 
further  simplified  by  a  purely  linear  resistive  network  for 

*’  Present  address.  Honeywell  Systems  and  Research  Center,  3640  Technol¬ 
ogy  Drive,  Minneapolis,  MN  55418. 


excitation  frequencies  below  a  few  megahertz  and  by  a  pure¬ 
ly  linear  capacitive  network  for  higher  frequencies.’’ 

The  purpose  of  the  present  paper  is  to  demonstrate  that 
an  equivalent  circuit  model  of  a  plasma  reactor  can  be  devel¬ 
oped  using  known  physical  properties  of  plasma  sheaths. 
The  main  feature  of  the  proposed  model  is  the  inclusion  of 
the  nonlinear  dep>endence  of  the  sheath  properties  on  the 
voltage  across  the  sheath.  The  purpose  of  the  model  is  to 
evaluate  the  effects  of  the  nonlinear  sheath  properties  on  the 
resulting  voltage  waveforms  and  the  dc  voltage  biases  across 
the  plasma  sheaths.  The  resulting  sheath  voltage  waveforms 
can  then  be  used  to  more  accurately  predict  the  ion  bom¬ 
bardment  energy  distribution  at  the  surface  of  the  electrodes 
in  a  planar  rf  plasma  reactor. 

II.  MODEL  OF  THE  PLASMA  SHEATH  IN  A  PLANAR  rf 
GASEOUS  DISCHARGE 

The  schematic  representation  of  a  planar  plasma  sheath 
used  in  the  following  discussion  is  shown  in  Fig.  2.  During 
the  discussion  of  the  sheath  properties  the  Bohm ' '  criterium 
will  be  used.  This  criterion  implies  that,  for  a  surface  biased 
strongly  negative  with  respect  to  a  plasma,  the  ions  are  accel¬ 
erated  across  a  pre-sheath  by  the  applied  electric  field  such 
that  they  enter  the  plasma  sheath  with  a  velocity  p>erpendicu- 
lar  to  the  sheath  boundary  and  with  a  magnitude  equal  to  or 
greater  than  Here,  k  is  Boltzmann's  constant, 

M  is  the  ion  mass,  and  7*,  is  the  electron  temperature. 

Several  assumptions  will  be  made  in  order  to  simplify 
the  present  discussion.  Refinements  of  the  model  can  be  in- 
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FIG.  I.  Schematic  representation  of  a  planar  rf  reactor. 
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FIG  2  Schematic  representation  of  a  planar  plasma  sheath 


eluded  in  future  discussions.  The  assumptions  used  for  the 
present  model  are  the  following: 

( 1 )  The  charged  particles  are  electrons  and  one  type  of 
positive  ions.  The  presence  of  more  than  one  type  of  positive 
ion  will  only  slightly  modify  the  model  presented.  In  many 
applications  of  plasma  processing  of  materials  electronega¬ 
tive  gases  are  used.  Plasmas  produced  in  these  gases  will 
contain  negative  ions.  The  presence  of  negative  ions  is  ex¬ 
cluded  in  the  present  model  since  an  accurate  incorporation 
of  these  ions  complicates  the  model  of  the  plasma  sheath. 
However,  as  long  as  the  negative  ion  density  inside  the  plas¬ 
ma  sheath  is  small  compared  to  that  of  the  positive  ion  den¬ 
sity,  the  present  model  is  a  good  approximation. 

(2)  The  electrons  and  ions  have  Maxwellian  velocity 
distributions  inside  the  plasma  with  temperatures  of  T,  and 
Ti ,  respectively. 

(3)  The  electron  temperature  T,  is  constant  in  time, 
while  the  ion  temperature  T,  is  equal  to  the  constant  neutral 
gas  temperature  T, .  The  assumption  of  constant  T,  limits 
the  model  to  a  minimum  frequency  of  the  applied  power 
source.  The  numerical  calculations  presented  in  this  paper 
for  a  planar  plasma  reactor  are  for  a  plasma  produced  in 
argon.  For  these  conditions,  the  energy  modulation  of  the 
majority  of  the  electrons  is  small  for  frequencies  larger  than 
approximately  SO  kHz  and  gas  pressures  smaller  than  ap¬ 
proximately  1  Torr.'^ 

(4)  The  electron  and  ion  densities  inside  the  plasma 
volume  are  constant  in  time.  The  validity  of  this  assumption 
depends  on  the  rate  of  loss  of  the  charged  particles  from  the 
plasma  volume  as  well  as  on  the  mechanism  by  which  the 
plasma  is  maintained.  The  charged  particle  density  modula¬ 
tion  due  to  loss  by  ambipolar  diffusion  in  a  plasma  reactor  is 
small  for  frequencies  larger  than  approximately  SO  kHz  and 
pressures  larger  than  approximately  0. 1  Torr  The  periodic 
charged  particle  current  across  the  plasma  sheaths  towards 
the  electrodes  is  governed  by  the  properties  of  the  plasma 
sheath.  The  model  shows  that  this  number  loss  of  charged 
particles  is  small  compared  to  the  charge  particle  density 
inside  the  plasma,  since  only  a  relatively  small  fraction  of  the 
plasma  electrons  and  ions  is  able  to  reach  the  electrodes. 

The  modulation  in  the  rate  of  production  of  charged 
particles  depends  on  the  maintenance  mechanism  of  the  dis¬ 
charge  plasma.  For  instance,  if  the  plasma  electrons  (ions) 
are  mainly  produced  by  secondary  electrons  emitted  due  to 
ion  bombardment  of  the  electrodes,  the  rate  of  charged  parti- 
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cie  production  inside  the  plasma  can  be  considerably  modu¬ 
lated.  This  is  caused  by  the  energy  and  flux  modulation  of 
these  electrons  due  to  the  time-dependence  voltage  differ¬ 
ence  across  the  plasma  sheath.  The  possible  effect  of  this 
phenomenon  is  not  included  in  the  present  model. 

( 5 )  The  influence  of  collisions  between  charged  parti¬ 
cles  and  neutral  particles  inside  the  plasma  sheath  on  the 
charged  particle  motion  can  be  neglected  (collisionless 
sheath).  For  instance,  this  is  a  good  assumption  for  the 
charged  particle  densities  (sheath  dimensions)  and  gas  pres¬ 
sures  (charged  particle  mean  free  paths)  used  for  anisotrop¬ 
ic  etching  of  rilms  in  a  planar  rf  plasma  reactor. 

(6)  The  transit  time  r,  of  ions  through  the  plasma 
sheath  is  not  larger  than  1//,  where /is  the  frequency  of  the 
applied  voltage.  If  this  condition  is  satisfied,  then  the  thick¬ 
ness  of  the  plasma  sheath  will,  at  each  time,  correspond  to 
the  steady-state  thickness  related  to  the  instantaneous  vol¬ 
tage  y^{t)  across  the  sheath.  Calculations  for  argon  ions  and 
plasma  densities  of  10‘°  cm“^  show  that  r,  is  smaller  than 
0.5  X  10“‘  s.  Therefore,  this  assumption  appears  to  be  valid 
for  frequencies  up  to  approximately  1  MHz. 

( 7 )  The  flux  of  electrons  collected  by  the  electrodes  de¬ 
pends  only  on  the  instantaneous  voltage  across  the  plasma 
sheath  and  not  on  the  frequency  of  the  voltage  waveform. 
Numerical  calculations  by  the  authors  indicate  that  this  as¬ 
sumption  is  valid  for  frequencies  up  to  a  few  megahertz. 

( 8 )  The  Bohm  sheath  criterium  can  be  used  in  the  form 
that  the  speed  with  which  the  ions  enter  the  plasma  sheath  is 
equal  to  (kT,/My’^  and  is  valid  for  the  nnge  of  voltages 
across  the  plasma  sheaths  of  interest.  ” 

(9)  The  electric  field  is  zero  at  the  boundary  between 
the  pre-sheath  and  the  plasma  sheath.'^ 

(10)  The  electrons  have  a  Boltzmann  density  distribu¬ 
tion  in  the  pre-sheath  and  the  plasma  sheath.'^ 

The  two  quantities  related  to  the  plasma  sheath  which 
have  to  be  included  in  the  equivalent  circuit  model  of  a 
planar  rf  plasma  reactor  presented  in  the  next  section  are  the 
capacitance  C,  of  the  plasma  sheath  and  the  conduction  cur¬ 
rent  /  through  the  sheath.  An  expression  for  C,  as  a  function 
of  sheath  potential  ( F,  —  F, )  can  be  derived  by  noting  that 
the  displacement  current  Ij  for  the  sheath  capacitance  is 
given  by'* 


r  _dQ_  dQ  dv,  _  ^  dV. 

“  dt  dV,  dt  '  dt  ' 


(1) 


where  Q  is  the  surface  charge  on  the  electrode.  For  a  planar 
electrode,  the  sheath  capacitance  C,  is  thus  given  by 


’  dA 


(2) 


where  fo  fhe  permittivity  of  free  space  and  it  is  assumed 
that  the  electric  field  at  the  electrode  surface  is  normal  to 
and  uniform  over  the  electrode  surface  with  area  /4. 

Using  the  mentioned  assumptions,  it  can  be  shown  that 


i/j 
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for  -  00  <  <  K,  (3a) 

£,<=0  for  (3b) 

Here, 

_  (A7',/2e)  (4) 

.u 

n,  =  n(,exp(ey,/kr,)  ,  (5) 


where  ng  is  the  density  of  electrons  and  ions  inside  the  plas¬ 
ma.  Equations  (2)  and  (3)  show  that  C,{y,)  is  time  vary¬ 
ing  and  has  a  strong  nonlinear  dependence  on  y,(l). 

The  assumptions  made  imply  that  the  conduction  cur¬ 
rent  /( y, )  through  the  plasma  sheath  is  given  by  the  static 
plasma  probe  current-voltage  characteristic  (Langmuir 
probe  curve)  derived  by  including  the  Bohm  criterium,  i.e., 

/(  K. )  =  AJ,  exp(ey,/kT. )  -b  AJ,  for  y,  (/,)<0. 

(6) 

Here,  A  is  the  area  of  the  electrode  and  the  saturated  charged 
particle  current  densities  are 

J,  =  -en„u,/4=  ~  eng(kT,/liTm)''^  (7) 

and 

J,  (8) 

where  m  and  M  are  the  electron  and  ion  mass,  respectively, 
and  V,  is  the  average  electron  speed  inside  the  plasma. 

The  expressions  given  in  this  section  for  the  capacitance 
C,  ( F, )  of  the  plasma  sheath  and  the  conduction  current 
/( y, )  through  the  plasma  sheath  will  be  used  in  the  equiva¬ 
lent  circuit  model  for  a  planar  rf  discharge  presented  in  the 
following  section. 

III.  EQUIVALENT  CIRCUIT  MODEL  FOR  A  PLANAR  rf 
PLASMA  REACTOR 

A  schematic  of  the  equivalent  electric  circuit  model 
used  in  the  present  discussion  is  shown  in  Fig.  3.  Here,  is 
the  voltage  of  the  applied  rf  signal  from  a  matched  rf  power 
supply;  y^  and  are  the  potentials  of  the  target  electrode 
and  the  plasma,  respectively,  and  Fss  =  and 
are  the  voltages  across  the  substrate-plas¬ 
ma  sheath  and  the  target-plasma  sheath,  respectively.  The 
blocking  capacitor  is  represented  by  ;  Csr  and  Jr  repre¬ 
sent  the  capacitance  of  and  conduction  current  through  the 


FIG.  3  Equivtlenl  electric  circuit  model  of  i  plasma  reactor. 
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sheath  adjacent  to  the  target  electrode,  respectively,  while 
C ss  and  Is  represent  the  corresponding  values  for  the  sheath 
adjacent  to  the  substrate  electrode.  The  expressions  for  Cgr 
and  C  ss  are  given  by  Eqs.  ( 2 )  and  ( 3 ) ,  while  the  expressions 
for/r  and  /j  are  given  by  Eqs.  (6),  (7),  and  (8). 

The  electrical  resistance  of  the  plasma  has  been  neglect¬ 
ed  in  this  circuit  model.  This  resistance  is  small  with  respect 
to  the  sheath  resistance  for  the  plasma  electron  densities  and 
voltage  frequency  range  considered  in  the  present  discus¬ 
sion.  However,  inclusion  of  the  plasma  resistance  does  not 
introduce  any  complications  for  the  circuit  model. 

The  purpose  of  the  model  is  to  derive  the  time  depen- 
denceofFp,  F^-,  the  total  current  at  the  node  denoted  by  F^^, 
and  the  conduction  currents  through  the  plasma  sheaths  in 
this  circuit.  From  these  results,  quantities  such  as  the  sheath 
potentials  as  a  function  of  time  and  the  resulting  dc  self  bi¬ 
ases  can  also  be  calculated. 

From  this  equivalent  circuit  model  and  from  conserva¬ 
tion  of  current  at  the  nodes  denoted  by  F^  and  F^ ,  it  follows 
that 

CsT  )  +  Cs  -^(  Frf  -  F, )  +  y,  =  0  ,  (9a) 

dt  dt 

Css  4-^p  +  Csr  )  +  /r  +  4  =  0  .  (9b) 

dt  dt 

If  Eqs.  (2)  and  (3)  for  C5(  F^ )  and  Eq.  (6)  for /(  F^ ) 
are  used  in  the  current  continuity  Eqs.  ( 9 ) ,  the  set  of  contin¬ 
uity  equations  can  be  solved  numerically  for  F^(f)  and 
Fy  ( / ) .  In  addition,  by  then  using  these  results,  the  current  in 
the  circuit  at  the  node  denoted  by  F,  and  the  conduction 
currents  /.v  and  7/  in  the  plasma  sheaths  can  be  calculated. 

IV.  NUMERICAL  RESULTS  AND  DISCUSSION 

Numerical  simulations  were  performed  based  on  the 
equivalent  circuit  model  for  the  following  assumed  condi¬ 
tions: 

M=  40  amu  (argon  discharge), 

=  10'”  particles/cm’  (charged  particle  density  in 
plasma), 

Ti  =  500  K  (ion  temperature), 

r,  =  23  200  K  (electron  temperature), 

Ar  =  100  TT  cm^  (area  of  target  electrode), 

Frf  =  10(X)  sin  (ot  (applied  rf  voltage). 

These  conditions  imply  that 

J,  =  214^A/cm^  (saturated  ion  current), 

y,  r=  —  38  mA/cm^  (saturated  electron  current), 

Ff  =  —  10.35  V  (floating  potential). 

The  floating  potential  is  the  voltage  difference  across  the 
plasma  sheath  and  pre-sheath  for  which  the  total  conduction 
current  through  the  sheath  is  zero.  Its  value  is  given  by 

Fy  =  -  {kT,/e)  ln(  -  JJJ, ) .  (10) 

which  follows  directly  from  Eq.  (6). 

Figure  4  shows  the  calculated  waveforms  for  the  voltage 
Fy  across  the  plasma  reactor  and  the  plasma  potential  F^ . 
'The  calculations  relate  to  a  frequency /=  1(X)  kHz,  a  block¬ 
ing  capacitance  Cg  =  150  pF,  and  an  electrode  area  ratio 
Ar/A,  =  1.0,  where  .4y  is  the  area  of  the  target  electrode 
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FIG.  4.  Calcuiated  waveforms  of  the  voltage  scross  the  plasma  reactor 
and  the  plasma  potential  V,  for  equal  areas  of  the  target  and  substrate  elec¬ 
trodes. 

and  A,  is  the  area  of  the  substrate  electrode.  The  sheath 
voltage  waveforms  f'ss  ^st  = 

shown  in  Fig.  S.  These  calculated  waveforms  show  that  both 
electrodes  always  have  a  negative  potential  with  respect  to 
the  plasma.  These  results  are  in  very  good  agreement  with 
measurements  by  Bruce.  ’’  Figure  6  shows  an  expanded  view 
of  the  waveforms  of  the  sheath  voltages  close  to  the  floating 


FIG.  5.  Csiculated  wiveform  of  the  poientUI  across  the  plasma  sheath 
adjacent  to  the  target  electrode  and  that  of  the  potential  y„  across  the  sub¬ 
strate  electrode  sheath  for  ^4  7. //4j  =  I. 


FIG.  6.  Calculated  sheath  voluge  waveforms  Fjy  and  Fjj  close  to  the  float¬ 
ing  potential  (OTAr/Ag  =  1. 

potential  Vp  where  it  is  seen  that  the  smallest  magnitude  of 
the  time  dependent  sheath  potential  is  only  slightly  smaller 
than  \  Vp\.  The  nonsinusoidal  nature  of  the  sheath  voltage 
waveforms  V and  V st  ( and,  hence,  of  K,  and  F,- )  is  due  to 
the  highly  nonlinear  properties  of  the  sheath  capacitances 
Css  and  C ST  and  the  conduction  currents and/^-  through 
the  plasma  sheaths.  At  this  frequency,  these  waveform 
shapes  are  predominantly  determined  by  the  large  difference 
in  magnitude  between  the  saturated  ion  current  and  the  larg¬ 
er  saturated  electron  current  and  by  the  exponential  increase 
of  the  electron  current  through  the  plasma  sheath  with  de¬ 
creasing  values  of  I  Fs  |.  Since  no  net  dc  current  component 
can  flow  in  the  circuit,  the  magnitude  of  the  sheath  poten¬ 
tials  need  only  be  slightly  smaller  than  |  Vp  |  for  part  of  the  rf 
cycle,  during  which  time  sufficient  net  electron  current  flows 
in  order  to  balance  the  net  ion  current  which  flows  during 
the  remainder  of  the  rf  cycle. 

The  waveform  of  the  conduction  current  Ip  through  the 
sheath  adjacent  to  the  target  electrode  is  shown  in  Fig.  7. 
This  current  becomes  independent  of  the  sheath  voltage  F 
(and  thus  time)  for  large  negative  sheath  voltages  (satu¬ 
rated  ion  current  only),  while  the  contribution  of  the  elec¬ 
tron  current  only  becomes  important  for  values  of  F^t  close 
to  the  floating  potential  Vp  (compare  with  Figs.  5  and  6). 
For  this  symmetrical  condition  (Ap  =  ,<5 )  tit®  waveform  of 
Isis  identical  with  that  of  Ip,  but  is  shifted  by  ir/o)  s,  as  can 
be  expected.  The  waveform  of  the  total  current,  /,o„i, 
through  the  discharge  is  shown  in  Fig.  8.  This  waveform  is 
again  nonsinusoidal  as  can  be  expected  from  the  nonlinear 
properties  of  the  plasma  sheaths.  It  should  also  be  noted  that 
the  time  averages  of  the  currents  Ip,  Is,  and  /,o„,  are  equal  to 
zero,  as  is  required,  since  no  net  dc  current  component  can 
flow  in  the  circuit. 
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TIME  (sec) 

FIG.  7.  Calculated  waveform  of  the  conduction  current  Ij-  through  the 
plasma  sheath  adjacent  to  the  target  electrode  for  A  j- /As  =  I. 


The  next  sequence  of  figures  gives  the  various  wave¬ 
forms  for  an  experimental  condition  in  which  the  area  of  the 
substate  electrode  is  five  times  the  area  of  the  target  electrode 
(Ar  =  IOOjt  cm^  and  =  SOOrr  cm^).  The  value  of  the 
blocking  capacitor  was  slightly  changed  in  these  calculations 
in  order  to  keep  the  voltage  amplitude  across  the  reactor, 
approximately  the  same  as  in  the  previous  calculations 
for  the  symmetrical  condition  (i4  7- //4s  =  1).  Figure  9  shows 


TIME  (Mc) 


FIG.  8.  Calculated  waveform  of  the  total  current  through  the  piaama 
reactor  for /4 7 //fj  —  I. 


TIME  (tec) 


FIG.  9.  Calculated  waveform  of  the  potential  across  the  plasma  sheath 
adjacent  to  the  target  electrode  and  that  of  the  potential  across  the  sub¬ 
strate  electrode  sheath  for  Aj-/As  =  0.2. 

the  waveforms  of  the  sheath  potentials  Vgr  and  Kjs ,  while 
Fig.  10  shows  an  expanded  view  of  these  waveforms  close  to 
the  floating  potential  The  maximum  value  of  |  Vgr  \  is 
about  twice  the  amplitude  of  the  voltage,  Vr,  across  the  plas¬ 
ma  reactor,  while  the  value  of  remains  close  to  the  value 
of  Vp.  The  latter  effect  is  shown  in  more  detail  in  Fig.  11. 


TIME  (sac) 


FIO.  10.  Calculated  sheath  voltage  waveforms  y„  and  close  to  the 
floating  potential  F^  fotAr/As  —  0.2. 
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FIG.  1 1.  Calculated  sheath  voltage  waveform  F $$  very  close  to  the  floating 
potential  Vr  fotAj-ZAs  =0.J. 

The  waveforms  of  the  conduction  currents,  Ij-  and  Is, 
through  the  two  plasma  sheaths  are  shown  in  Figs.  12  and 
13,  respectively.  The  waveform  of/j-  is  analogous  to  that  for 
the  symmetrical  configuration.  However,  the  electron  cur¬ 
rent  contributes  significantly  to  Is  during  all  times,  since  the 
saturated  ion  current  at  the  substrate  electrode  is  336  mA, 
while  at  the  target  electrode  this  current  is  only  67  mA.  The 


FIG.  12.  Calculated  waveform  of  the  conduction  current  It  through  the 
plasma  sheath  adjacent  to  the  target  electrode  fot  At/As  0.2. 


FIG.  1 3.  Calculated  waveform  of  the  conduction  current  through  the 
plasma  sheath  adjacent  to  the  substrate  electrode  fotAr/A,  »  0.2. 

reason  that  the  amplitude  of  » *"“ch  larger  than  that  of 
Fss  is  the  difference  in  electron  and  saturated  ion  currents  at 
the  two  electrodes  due  to  the  difference  inAr  and  .,4s .  The 
variation  in  Fjs  required  in  order  to  adjust  Is  for  the  vari¬ 
ation  in  Ij-  becomes  smaller  than  the  variation  in  Kjt  the 
ratio  Aj-ZAs  becomes  smaller.  The  waveform  of  the  total 
current  through  the  plasma  reactor  is  shown  in  Fig.  14. 

V.  CONCLUSIONS 

The  model  presented  is  based  on  the  inclusion  of  the 
physical  properties  of  plasma  sheaths  in  an  equivalent  circuit 
representation  of  a  planar  rf-excited  gaseous  discharge.  The 


FIG.  14.  Calculated  waveform  of  the  total  current  through  the  plaaroa 

reactor  tor  At/A j  =  0.2. 
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proposed  model  is  valid  under  the  assumptions  made  in  Sec. 
II.  The  model  shows  that  the  waveforms  of  the  voltage  dif¬ 
ferences  across  the  plasma  sheaths  are  highly  nonsinusoidal. 
These  waveforms  are  in  agreement  with  studies  reported  by 
Bruce.”  The  resulting  current  waveforms  in  the  circuit  arc 
also  highly  nonsinusoidal.  The  waveforms  of  the  voltage, 
Fr  •  across  the  reactor  and  of  the  total  current,  /  ,o,,, ,  through 
the  reactor  have  important  implications  when  determining 
the  power  applied  to  the  reactor.  This  model  also  shows  that 
the  relative  amplitudes  of  the  voltage  waveforms  across  the 
target  and  substrate  plasma  sheaths  depend  on  the  area  ratio 
of  the  electrodes,  with  the  largest  amplitude  waveform 
across  the  sheath  adjacent  to  the  smallest  area  electrode. 

The  effect  of  the  calculated  nonsinusoidal  waveform  of 
the  plasma  sheath  voltage  on  the  energy  distribution  of  the 
positive  ions  arriving  at  the  electrodes  has  been  investigated 
and  will  be  reported  in  a  separate  paper.  An  extension  of  the 
model  to  include  the  presence  of  negative  ions  inside  the 
plasma  sheath  is  currently  in  progress. 
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'The  results  described  in  this  paper  are  also  relevant  for  plasma  reactors  of 
different  geometries  as  long  as  the  thickness  of  the  plasma  sheaths  is  smalt 
compared  to  the  electrode  dimensions.  Calculations  show  that  for  the  elec¬ 
tron  (ion)  number  densities  typical  in  low-pressure  plasma  reactors,  the 
thickness  of  the  plasma  sheath  is  of  the  order  of  1  mm.  The  model  calcula¬ 
tions  presented  here  can  also  be  extended  to  take  into  account  the  effect  of 
reactor  walls,  etc.,  on  the  electrical  characteristics  of  the  reactor  by  adding 
circuit  elements  to  the  equivalent  electric  circuit  model  discussed  in  Sec. 
III. 
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ABSTRACT 

A  model  Is  presented  for  the  energy  distribution  of  ions  bombarding 
electrode  surfaces  in  planar  rf  plasma  reactors  at  low  pressures,  vdiere 
collisions  in  the  sheath  may  be  ignored.  Given  the  waveform  for  the  voltage 
difference  across  the  sheath  adjacent  to  an  electrode,  the  model  predicts 
the  ion  bombardment  energy  distribution  in  the  limit  ti«l/f  and  tr»l/f, 
where  tr  is  the  transit  time  of  the  ions  across  the  sheath  and  f  is  the 
frequency  of  the  applied  rf  voltage.  Utilizing  sheath  voltage  waveforms  from 
a  previously  published  equivalent  circuit  model  of  a  rf  reactor^,  the  model 
shows  that  the  ion  bombardment  energy  distribution  is  blmodal  and  skewed 
toward  lower  ion  bombardment  energies.  The  model  also  demonstrates  one 
of  the  mechanisms  responsible  for  narrowing  of  the  ion  bombardment 
energy  distribution  with  increasing  rf  frequency.  These  results  are 
compared  with  experimental  measurements. 

I.  INTRODUCTION 

Radio  frequency  gas  discharges  currently  play  a  major  role  in  the 
microelectronics  industry  for  the  fabrication  of  integrated  circuits.  These 
discharges  are  typically  used  for  stripping  of  photoresist,  anisotropic 
etching  of  semiconductor,  oxide,  and  metal  surfaces,  and  the  deposition  of 
organic  and  Inorganic  thin  films.  In  a  rf  discharge,  positive  ions  produced  in 
the  plasma  volume  are  accelerated  across  the  plasma  sheath  and  may  arrive 
at  the  electrodes  with  significant  bombardment  energy.  Because  the  sheath 


‘^Present  address:  Honeywell  Systems  and  Research  Center.  3660  Technology  Drive. 
Minneapolis.  MN  55418. 

Metze,  D.  W.  Emle.  and  H.  J.  Oskam.  J.  Appl.  Phys.  60.  3061  (1986). 
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structure  (e.g..  sheath  thickness,  sheath  voltage,  etc.)  In  an  rf  discharge  is  a 
periodic  function  of  time,  the  resulting  time  averaged  bombardment  energy 
is  not  constant,  but  has  a  distribution  of  bombardment  energies.  In  all  the 
above  mentioned  applications  .  the  ion  bombardment  energy  distribution 
with  which  ions  bombard  the  surfaces  being  processed  is  a  critical  factor  in 
determining  the  resulting  device  quality  and  properties.  For  example,  the 
ion  bombardment  energy  distribution  determines  the  degree  of  anisotropy 
in  thin  film  etching,  the  amoimt  of  ion  impact  induced  surface  d’unage  to 
semiconductor  and  oxide  surfaces,  the  erosion  rate  of  photoresist,  and  the 
degree  of  cross  linking  In  polymer  film  deposition.2.3  Therefore,  it  is 
important  to  obtain  a  better  understanding  of  the  dependence  of  the  ion 
bombardment  energy  distribution  on  discharge  parameters  such  as  reactor 
geometry,  frequency  of  the  applied  rf  voltage,  gas  pressure,  gas  composition, 
etc. 

The  ion  bombardment  energy  distribution  in  rf  plasma  reactors  has  been 
investigated  experimentally  by  several  authors.  Vasile  and  Smolinsky 
studied  the  ion  bombardment  energy  distribution  at  the  wall  of  a  13.6  Mhz  rf 
discharge  in  a  mixture  of  argon  and  vlnyltrlmethylsilane  at  a  few  hundred 
milllTorr  utilizing  a  retarding  grid  energy  analyzer  followed  by  a  quadrupole 
mass  spectrometer.^  A  similar  arrangement  was  used  by  Thompson  et  al  to 
measure  the  ion  bombardment  energy  distribution  at  an  electrode  of  a  13.6 
MHz  discharge  in  SFe,  CF3CI,  and  CFaBr  for  pressures  of  0.2- 1 .0  Torr.®  In 
both  cases,  the  measured  ion  bombardment  energy  distribution  consisted  of 
a  single  peak  at  a  bombardment  energy  and  with  an  energy  spread 
determined  by  the  discharge  pressure  and  power.  Measurements  of  the  ion 
bombardment  energy  distribution  at  the  electrodes  of  a  planar  rf  reactor 
have  been  studied  by  Cobum  and  colleagues  for  frequencies  of  100  KHz- 13.6 
MHz  in  Ar  at  pressures  of  20-50  mTorr  using  an  electrostatic  energy 


^D.  L.  Flamm,  V.  M.  Donnelly,  and  D.  E.  Ibbotson,  in  VLSI  Electronics,  edited  by  N.  G. 
Einspruch  and  D.  M.  Brown  (Academic.  Orlando.  1984).  VoL  8.  Chap.  8.  p.  189. 

^Russ.  A.  Moigan.  Plasma  Etching  In  Semiconductor  Fabrication.  Plasma  Technology.  1 
(Elsevier.  New  York,  1985). 

J.  Vasile  and  G.  Smolinsky.  Int.  J.  Mass  Spectrum.  Ion  Phys.  12.  133  (1973). 

E.  Thompson.  K.  D.  Allen.  A.  D.  Richards,  and  H.  H.  Sawln.  J.  ^pl.  Phys.  89.  1890  (1986). 
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anal3^er  followed  by  a  quadrupole  mass  spectrometer.6-®  Restilts  from  these 
studies,  in  which  collision  effects  were  minimized  by  a  Judicious  choice  of 
the  analyzed  ion  species,  showed  that  the  shape  of  the  ion  bombardment 
energy  distributions  is  a  function  of  the  rf  frequency,  the  ion  mass,  and  the 
gas  pressure.  The  shape  was  bimodal  for  low  gas  pressure,  low  rf  frequency, 
and  low  mass  ions,  but  became  a  single  peak  at  higher  pressure,  frequency, 
and  ion  mass.  Similar  results  have  recendy  been  reported  by  Brland  et  oL  in 
Ar,  CF4,  Eind  SFe  rf  discharges  9 

Theoretical  work  has  also  been  performed  in  an  effort  to  calculate  the  ion 
bombardment  energy  distribution  in  a  planar  rf  reactor.  Tsui  modeled  the 
ion  bombardment  energy  distribution  in  a  rf  reactor,  calculating  a  bimodal 
distribution  shape,  with  the  distribution  skewed  toward  higher  energies. 
Okamoto  and  Teimagawa  subsequendy  calculated  the  energy  spread  of  the 
ion  bombardment  energy  distribution  at  the  electrode  of  an  rf  reactor  as  a 
function  of  rf  frequency  in  the  limit  where  the  transit  time  of  the  ions 
through  the  sheath,  tr.  is  large  with  respect  to  the  period  of  the  applied  rf 
voltage,  1/f,  where  f  is  the  frequency  of  the  applied  voltage.^!  Their  results 
indicated  that  the  ion  bombardment  energy  distribution  is  single  peaked. 
Recendy,  Suzuki  et  al.  have  have  considered  the  limit  where  tr«l/f.^2  Th^ 
calculated  the  minimum  and  maximum  bombardment  energy  and  the 
spread  of  the  energy  distribution.  While  all  of  the  above  authors  ignored  the 
effect  of  collisions  in  the  sheath,  Monte  Carlo  techniques  have  recendy  been 
used  by  Kushneri®  and  Thompson^^  et  al.  to  investigate  the  ion 
bombeirdment  energy  distribution  at  the  electrode  of  a  planar  rf  plasma 
reactor  in  the  presence  of  collisions.  The  Monte  Carlo  simulations  show  the 

®J.  W.  Cobum  and  E.  Kay,  J.  Appl.  Phys.  43,  4965  (1972). 

^K.  KOhler,  J.  W.  Cobum.  D.  E.  Home,  and  E.  Kay,  J.  ^pL  PltyS.  57,  59  (1985). 

®K.  KChler,  D.  E.  Home,  and  J.  W.  Cobum,  J.  Appl.  Pbys.  58, 3350  (1985). 

9p.  Brland,  G.  Turban,  and  B.  Grollean,  Mat.  Res,  Soc.  Symp.  Proc.  68.  109  (1986). 

IOrt.C.  Tsui,  Phys.  Rev.  168.  107  (1968). 

i  iy.  Okamoto  and  H.  TamagaM^,  J.  Pl^.  Soc.  Japan  29. 187  (1970). 

^^K.  Suzuki,  K.  Nlnomiya,  S.  Nlshknatsu,  and  S.  Okudalra,  J.  Vac.  Scl.  Technol.  B3.  1025 
(1985). 

J.  Kushner,  J.  Appl.  Phys.  58,  4024  (1985). 

E.  Thompson.  H,  H.  Sawln,  and  D.  A.  Fisher,  (private  communications). 


dependence  of  the  shape  of  the  distribution  on  discharge  parameters  such 
as  rf  frequency,  ion  mass,  pressure,  etc.  The  results  display  both  single 
peak  and  bimodal  ion  bombardment  energy  distributions,  depending  on  the 
discharge  conditions.  All  of  the  above  studies  were  performed  assuming  that 
the  electric  potential  difference  across  the  plasma  sheath  is  the 
superposition  of  a  dc  bias  voltage  and  a  sinusoidal  rf  voltage,  or  a  version  of 
this  shape  clipped  to  the  floating  potential  of  the  sheath. 

This  paper  presents  calculations  of  the  ion  bombardment  energy 
distribution  at  the  electrode  of  a  planar  rf  plasma  reactor.  A  theoretical 
model  is  developed  in  the  next  section  for  the  ion  bombardment  energy 
distribution  at  a  surface  exposed  to  the  plasma.  The  model  assumes  a 
collisionless  plasma  sheath.  The  application  of  this  model  in  the  limits 
tr«l/f  and  tr»l/f  will  be  discussed.  Then,  utilizing  sheath  voltage 
waveforms  predicted  by  an  equivalent  circuit  model  for  a  planar  rf  plasma 
reactor.  1  calculations  will  be  presented  for  the  ion  bombardment  energy 
distribution  at  the  electrode  for  both  limiting  cases.  The  effect  of  the 
various  discharge  parameters  on  the  ion  bombardment  energy  distribution 
will  be  discussed  and  the  results  compared  with  experiments.  The 
resulting  ion  bombardment  energy  distribution  has  Important  consequences 
for  the  fabrication  of  microelectronic  devices. 


n.  MODEL  OF  THE  ION  BOMBARDMENT  ENERGY  DISTRIBUTION 


A.  Plasma  Sheath  Model 

The  model  for  a  planar  plasma  sheath  in  a  rf  discharge  used  in  this  paper 
is  similar  to  the  plasma  sheath  model  presented  previously.  ^  A  schematic 
representation  of  the  plasma  sheath  is  shown  in  Fig.  1 .  In  order  to  arrive  at 
a  tractable,  yet  physically  realistic  solution,  the  following  assumptions  will  be 
made.  For  a  detailed  discussion  of  the  applicability  of  these  assumptions 
refer  to  Ref.  1. 

(1)  The  plasma  consists  of  a  singly  charged  positive  ion  species, 
electrons,  and  neutrals. 

(2)  The  ions  and  electrons  have  a  Maxwellian  energy  distribution  with 
temperatures  Ti  and  Tg,  respectively. 
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SURFACE  SHEATH  PRE- SHEATH  PLASMA 


Fig.  1  Schematic  representation  of  a  planar  plasma  sheath. 
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(3)  The  ion  and  electron  temperatures  Ti  and  Te  are  Independent  of 
time  and  spatial  coordinates,  with  Ti  being  equal  to  the  neutral  gas 
temperature. 

(4)  The  ion  and  electron  densities  Inside  the  plasma  volume,  Uq  and  n^, 
respectively,  are  constant  in  time. 

(5)  The  plasma  sheath  is  collisionless. 

(6)  The  ions  enter  the  sheath  from  the  pre-sheath  with  a  speed  equal  to 
(kTe/M)i/2,  as  required  by  the  Bohm  sheath  crlterlum.^s  Here,  k  is 
Boltzman's  constant  and  M  is  the  ion  mass 

(7)  The  electric  field  ^(x)  is  zero  at  the  boundary  between  the  pre-sheath 
and  the  plasma  sheath. 

(8)  The  electrons  have  a  Boltzman  density  distribution  in  the  pre-sheath. 

(9)  The  influence  of  the  electrons  in  the  sheath  can  be  neglected.  The 
validity  of  this  assumption  depends  on  the  voltage  across  the  plasma  sheath 
and  the  sheath  propc'itles  of  Interest.  For  a  surface  biased  negatively  with 
respect  to  the  plasma  so  that  Vs«Vi,  vdiere  Vi=-  kTe/2e  and  Vs  is  the 
voltage  across  the  plasma  sheath,  only  a  small  fraction  of  the  electrons 
penetrate  a  significant  distance  into  the  sheath.  Thus  the  electrons  have 
little  effect  on  the  electric  field  in  the  majority  of  the  sheath  region.  As  a 
consequence,  the  electrons  have  little  effect  on  the  motion  of  the  ions 
through  the  plasma  sheath.  Since  this  is  the  parameter  of  interest,  as 
discussed  below,  the  influence  of  electrons  in  the  sheath  can  be  ignored  for 
situations  where  Vs«Vi.  Indeed,  this  is  always  the  case  for  an  electrode 
surface  in  a  planar  rf  reactor,  as  demonstrated  by  the  sheath  voltage 
waveforms  calculated  in  Ref.  1. 

(10)  The  rate  at  which  ions  enter  the  sheath  from  the  pre-sheath  is 
constant.  This  assumption  is  valid  in  the  limit  tr«l/f.  For  this  limit, 
although  the  sheath  thickness  d  will  be  a  function  of  time,  d  will  change 
slowly  relative  to  the  response  time  of  the  ions.  For  the  limit  tr»l/f,  this 
assumption  is  also  valid,  since  the  sheath  thickness  d  will  then  be  a  constant 
determined  by  the  time  averaged  sheath  voltage  Vdc  (l.e.,  the  dc  self  bias), 
since  the  response  time  of  the  ions  will  be  slow  relative  to  changes  in  the 
sheath  voltage.  Calculations  for  the  discharge  parameters  used  in  the 


Bohm,  in  The  Characteristics  of  Electrical  Discharges  In  Magnetic  Fields,  edited  by  A 
Guthrie  and  R  K,  Wakerllng  (McGraw-Hill,  New  York,  1949),  Chap.  3. 
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following  discussion  show  that  this  assumption  is  valid  for  frequencies  less 
than  approximately  1  MHz  and  greater  than  approximately  lOMHz. 

Using  the  above  assumptions.  Poisson's  equation  in  the  sheath  region  is 


d^2 

-  2eni 

f-  2eVi  v/2 

r-  2eV(x)>i/2 

dV(x)  ■ 

to 

1  M  J 

L  M  J 

where  the  density  of  ions  and  electrons  at  the  boundary  between  the  pre¬ 
sheath  and  the  sheath  is 


ni  =  no  exp(eVi/kTe). 

Under  assumption  (7),  Ekj.  (1)  can  be  integrated  to  yield 

^ivwi  =  -  2  ^  if  ^ 

Likewise,  since  4M=-  dV(x)/dx.  Eq.  (3)  can  be  integrated  to  3deld 
.  2(  to  >/2  ,  -1/2  r/V(x)Nl/2  -11/2  r/V(x)Nl/2  1 

Li-vrj  -iJ  [[vrj 

where  the  sheath  thickness  d  is  given  by 


(2) 


(3) 


(4) 


d  s  3 


2f  to  ^l/2 


(* 


(5) 


Using  E:qs.  (3).  (4).  and  (5).  it  can  be  shown  that 

I,,  ^  /^-eni,,  v/2  [Tb  b2  >il/211/3  fb  b2  xi/2-|i/3l 

Ux)  =  -2^— ViJ  |2J+(i+452)  J  J  }• 


(6a) 


where 


■iferc-v.) 


1/2 


and 


b  a  d  -  X 


(6b) 


(6c) 


The  above  results  for  Eqs.  (3) -(6)  are  equivalent  to  the  Child -Langmuir 
sheath  model,  where  Vi*-  kTe/2e,  as  required  by  the  Bohm  crlterlum.  The 
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above  expressions  give  d  and  ^(x)  as  a  function  of  the  voltage  Vs  across  the 
sheath. 


B.  Ion  Bombardment  Energy  Distribution  For  tr«l/f 

For  a  planar  rf  plasma  reactor,  the  sheath  voltage  is  a  function  of  time, 
Vs(t).  If  the  transit  time  of  the  ions  across  the  sheath  is  short  with  respect 
to  the  period  of  the  applied  rf  voltage,  l.e.,  tr«l/f,  the  sheath  parameters  d 
and  ^(x)  will  be  given  by  Eqs.  (5)  and  (6),  with  Vs=Vs(t).  Therefore,  d  and 
^(x)  are  also  functions  of  time,  l.e.,  d(t)  and  x(x,t).  As  a  consequence,  the 
resulting  ion  bombardment  energy  for  an  ion  E(t)  will  be  a  function  of  time 
and  will  be  determined  by  the  time  at  which  the  ion  enters  the  sheath 
region. 

Consider  the  following  equation  of  motion  for  an  ion  which  enters  the 
sheath  at  time  t=ti, 

M  ^x(t)  =  e  ^(x,t)  ,  (7a) 

with  initial  conditions 

x(t=ti)  =  d(t=ti)  (7b) 


and 


t=ti 


f-  2eViV/2 

I  M  J 


(7c) 


Equation  (7)  can  be  solved  for  the  velocity  of  the  ion  dx(t)/dt  at  any  time  t^q 
and  for  the  transit  time  of  the  ion  across  the  sheath,  tr(ti),  which  is  now  a 
function  of  the  time  when  the  ion  entered  the  sheath.  The  resulting 
bombardment  energy  E(t)  of  the  ion  is  then  given  by 


E(ti+tr(ti)l 


M 

2 


t=ti+tr(ti)  * 


(8a) 


However,  since  the  sheath  voltage  is  periodic,  all  the  sheath  parameters  and 
the  resulting  ion  bombardment  energy  will  be  periodic.  Therefore, 

E[ti+tr(ti)]  =  E(<|>)  .  (8b) 
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where  0  is  the  phase  angle  of  the  sheath  voltage  at  the  time  ti  and  where 

If  PE(E)dE  is  the  probability  that  an  ion  will  strike  the  surface  with  a 
bombardment  energy  In  the  Interval  dE  and  P0((]>)d<t>  Is  the  probability  that  an 
Ion  is  injected  Into  the  sheath  with  a  phase  angle  in  the  interval  d<t>,  then 

Pe(E)  =  P4,((|))  . 

As  a  result  of  assumption  (10).  P,^(4>)=l/2rc.  Therefore,  the  ion  bombardment 
energy  distribution  Pe(E)  Is  given  by 

pe® = ^  iHi  ■  O' 

where  |d())/dE|  Is  determined  from  Eq.  (8). 

It  is  worth  noting  that 


lim  E(<{»)  -  eVs«»)  , 

tff — >0 

where  Vs(4))  Is  the  sheath  voltage  at  phase  angle  <>.  Therefore,  Eq.  (9) 
simplifies  for  trf-^O  to 


(10) 


C.  Ion  Bombardment  Energy  Distribution  For  tr»l/f 

If  tr»l/f.  the  ions  will  not  be  able  to  respond  to  the  Instantaneous 
sheath  voltage,  but  only  to  the  average  value  V(jc.  In  this  situation,  the  sheath 
thickness  d  Is  constant  and  can  be  assumed  to  be  given  by  Eq.  (5).  with 
Vs=Vdc-  Likewise,  the  electric  field  in  the  sheath  can  be  assumed  to  have  a 
constant  shape,  but  a  time  dependent  magnitude,  as  given  by 


.  ,  ,  ,  Vs(t) 

£(x,t)  =  ^dc(x) 


(11) 
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where  ^dc(x)  Is  given  by  Eq.  (6)  with  Vs=Vdc.  For  this  case,  Eq.  (7),  (8).  and 
(9)  can  again  be  solved  utilizing  the  constant  sheath  thickness  assumption 
and  Eq.  (11),  yielding  the  ion  bombardment  energy  distribution  Pe(E). 


m.  NUBIERICAL  RESULTS  AND  DISCUSSION 

Numerical  calculations  were  performed  based  on  the  above  ion 
bombardment  energy  distribution  model  for  the  discharge  conditions 
assumed  in  Ref.  1,  i.e., 

M  =  40  amu  (argon  discharge), 

no  =  IQiO  partlcles/cm3  (charged  particle  density  in  the  plasma), 

Ti  =  500  OK  (ion  temperature), 

Te  =  23,200  oR  (electron  temperature). 

For  these  discharge  parameters,  calculations  show  that  0.1  ps  <  tr  <  0.5  ^is. 

The  sheath  voltage  waveforms  for  the  target  electrode  presented  in  Ref.  1 
were  used  for  Vs(t)  in  these  calculations.  These  waveforms  are  reproduced 
in  Figs.  2  and  3  for  future  reference.  Bg.  2  is  for  a  planar  rf  reactor  with  an 
electrode  area  ratio  At/As=1.0,  while  Fig.  3  is  for  At/As=0.2,  where  At  is 
the  area  of  the  target  electrode  and  As  is  the  area  of  the  substrate  electrode. 
These  sheath  voltage  waveforms  were  derived  for  f=100  KHz  and  can 
therefore  be  directly  used  in  the  ion  bombardment  energy  distribution 
model  for  the  limiting  case  tr«l/f.  Calculations  were  performed  at  f=100 
KHz  using  these  waveforms.  In  addition,  in  order  to  investigate  the  effect  of 
frequency  on  the  ion  bombardment  energy  distribution,  calculations  were 
also  performed  at  fsl3.6  MHz,  vdiere  the  limiting  situation  tr»l/f  applies. 
For  this  situation,  the  same  sheath  voltage  waveforms  from  Figs.  2  md  3 
were  again  used,  even  though  these  waveforms  probably  do  not  accurately 
reflect  the  actual  sheath  voltage  waveforms  at  this  frequency.  However,  as 
shown  below,  these  calculations  are  useful  for  demonstrating  the  effect  of 
frequency  on  the  resulting  ion  bombardment  energy  distribution, 
independent  of  the  effect  of  frequency  on  other  sheath  properties,  such  as 
the  sheath  voltage  waveform. 

In  order  to  more  accurately  compare  the  results  presented  below  with 
experimental  measurements  of  the  ion  bombardment  energy  distribution,  it 
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TIME  (sec) 


Fig.  2.  Calculated  waveform  of  the  potential  Vsr  across  the  plasma  sheath 
adjacent  to  the  target  electrode  (and  that  of  the  potential  Vss  across  the 
substrate  electrode)  for  At/As=1,  from  Ref.  1. 


CU  (U 


TIME  (sec) 

Fig.  3.  Calculated  waveform  of  the  potential  Vst  across  the  plasma  sheath 
adjacent  to  the  target  electrode  (and  that  of  the  potential  Vss  across  the 
substrate  electrode)  for  At/As=0.2,  from  Ref.  1 . 
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was  desirable  to  fold  in  the  effect  of  the  finite  energy  resolution  of  an 
experimental  energy  analysis  system.  If  it  is  assumed  that ,  to  first  order, 
the  energy  anatyzer  has  a  square  transmission  window  of  width  AE,  the 
quantity  measured  by  the  experimental  system  will  be  the  ion  bombardment 
energy  distribution  Pe(E)AE.  In  the  following  discussion,  results  are 
presented  for  Pe(E)AE  with  AE=0.5  eV. 

Figures  4  and  5  show  the  calculations  of  the  ion  bombardment  energy 
distribution  for  two  different  electrode  area  ratios  utilizing  the  target  sheath 
voltage  waveforms  from  Figs.  2  and  3.  respectively,  for  f=100  KHz  and  1 
MHz.  As  can  be  seen  from  these  results,  the  ion  bombardment  energy 
distributions  are  blmodal  and  skewed  toward  lower  energies.  The  blmodal 
characteristic  of  the  distributions  can  be  qualitatively  explained  by 
considering  Eq.  (10)  for  the  limit  trf->0.  In  this  limit.  Pe(E)AE  is  inversely 
related  to  the  derivative  with  respect  to  phase  angle  (or  time)  of  the  sheath 
voltage  waveform  .  Thus,  the  resulting  ion  bombardment  energy  distribution 
has  a  lower  value  for  those  energy  ranges  vdiich  correspond  to  locations  on 
the  sheath  voltage  waveform  vdiere  the  waveform  is  a  rapidly  changing 
function  of  phase  angle  (or  time).  In  addition,  the  skewing  of  the 
distributions  is  due  to  the  nonsinusoidal  shape  of  the  sheath  voltage 
waveforms  which  results  in  the  ions  being  accelerated  for  a  large  fraction  of 
each  rf  cycle  by  a  small  sheath  potential  and  only  for  a  small  fraction  of  each 
rf  cycle  by  a  large  sheath  potential.  Thus,  a  larger  number  of  ions  impact  the 
target  electrode  with  low  bombardment  energy  than  with  high 
bombardment  energy. 

Figs.  4  and  5  also  show  that  the  energy  spread  of  the  distribution  narrows 
with  increasing  frequency.  This  is  to  be  e3q}ected.  since  at  low  excitation 
frequencies  (e.g.,  at  100  KHz)  the  ions  traverse  the  sheath  in  a  fraction  of 
the  rf  period.  Consequently,  the  ions  bombard  the  the  sheath  with  an 
energy  -eVs(ti),  where  Vs(ti)  is  the  sheath  voltage  at  the  instant  when  the 
ions  enter  the  sheath.  Therefore,  the  resulting  energy  spread  in  the 
distribution  function  is  essentially  determined  by  the  peak  to  peak 
amplitude  of  the  sheath  voltage  waveform.  If,  however,  the  ion  transit  time 
is  long  with  respect  to  the  rf  period  (e.g.,  at  13.6  MHz),  then  the  ions  are 
effectively  accelerated  across  the  sheath  by  the  time  average  of  the  sheath 
voltage  waveform.  In  this  case,  a  smaller  percentage  of  the  energy  gained  by 
the  ions  in  the  sheath  can  be  attributed  to  the  rf  component  of  the  sheath 
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Fig.  4.  Calculated  ion  bombardment  energy  distribution  for  At/As=1. 


2 


Fig. 


ION  IMPACT  ENERGY  (eV) 


5.  Calculated  Ion  bombardment  energy  distribution  for  At/As=0.2. 
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voltage.  As  a  result,  the  energy  spread  of  the  Ion  bombardment  energy 
distribution  decreases  with  increasing  frequency. 

The  results  of  these  calculations  can  be  compared  with  the  experimental 
measurements  of  Cobum  and  Kay.^  In  their  studies,  the  energy  distribution 
of  the  ion  surface  flux  at  an  electrode  was  determined  for  low  pressures  £md 
for  minority  ion  species  vdilch  did  not  undergo  significant  resonance  charge 
transfer  collisions.  Thus,  the  experimental  conditions  closely  approximated 
the  collisionless  sheath  assumption  used  in  deriving  our  model.  The  shape 
of  the  energy  distributions  shown  in  Figs  4  and  5  agree  well  with  the 
experimental  results.  In  both  cases,  the  distributions  are  bimodal  and 
skewed  toward  lower  energy.  In  addition,  the  experimental  studies  showed 
the  distribution  narrowing  with  Increasing  ion  mass,  for  a  fixed  rf  frequency. 
This  narrowing  is  similar  to  the  narrowing  of  the  spread  in  the  distribution 
predicted  by  our  model.  This  similarity  results  from  the  fact  that  both 
effects  are  due  a  change  in  the  ratio  of  the  ion  transit  time  to  the  rf  period. 
This  ratio  may  be  increased  by  either  increasing  the  ion  mass,  which  causes 
an  increase  in  the  ion  transit  time,  or  by  increasing  the  rf  frequency,  which 
decreases  the  rf  period. 

Finally,  a  comparison  between  Figs.  4  and  5  shows  that  the  ion 
bombardment  energy  achieves  higher  values  in  Fig.  5  than  in  Fig.  4.  This 
occurs  since  the  area  of  the  target  electrode  is  smaller  than  the  area  of  the 
substrate  electrode  for  the  data  given  in  Fig.  5,  while  the  areas  are  equal  in 
for  Fig.  4.  Since  the  sheath  voltage  waveforms  in  Figs.  2  and  3  were  both  for 
the  same  peak  to  peak  rf  voltage  applied  across  the  reactor  electrodes,  the 
target  sheath  voltage  waveform  had  a  larger  peak  to  peak  magnitude  for  the 
case  of  unequal  area  electrodes.  Thus,  the  ion  bombardment  energies  eire 
higher  in  Fig.  5.  where  the  target  electrode  Is  smaller  than  the  substrate 
electrode. 


IV.  CONCLUSIONS 

A  model  has  been  presented  for  the  ion  bombardment  energy  distribution 
at  an  electrode  of  a  planar  rf  plasma  reactor.  The  model  is  valid  imder  the 
conditions  given  in  Sec.  II,  where  it  was  assumed  that  the  sheath  adjacent  to 
the  electrode  is  collisionless  and  that  the  electron  temperature  and  the 
charged  particle  density  in  the  plasma  are  independent  of  time.  Given  the 
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sheath  voltage  waveforms,  the  model  predicts  the  ion  bombardment  energy 
distribution  in  the  limit  tr«l/f  and  tr»l/f,  where  tr  is  the  transit  time  of 
the  ions  across  the  sheath  and  f  is  the  frequency  of  the  applied  rf  voltage. 
Numerical  calculations  were  performed  at  f^lOO  KHz  and  13.6  MHz  using 
this  model  with  the  sheath  voltage  waveforms  given  in  Ref.  1  for  the  target 
electrode.  Although  these  sheath  voltage  waveforms  are  strictly  valid  only 
for  the  calculations  at  f=100  KHz.  vdiere  tr«l/f.  they  were  also  used  for  the 
calculations  at  f=13.6  MHz,  where  tr»l/f.  in  order  to  determine  the  effect 
of  frequency  on  the  ion  bombardment  energy  distribution  independent  of 
the  effect  of  frequency  on  other  sheath  properties. 

The  resulting  numerical  calculations  show  that  the  ion  bombardment 
energy  distribution  is  blmodal  and  skewed  toward  lower  energy.  In 
addition,  the  energy  spread  of  the  distribution  narrows  with  increasing 
frequency.  Finally,  the  Increase  in  ion  bombardment  energy  with  decreasing 
target  electrode  to  substrate  electrode  area  ratio  is  demonstrated.  A 
qualitative  explanation  of  these  effects  has  been  given.  These  results 
compare  favorably  with  previous  experimental  studies.® 

As  Indicated,  the  above  results  were  derived  under  the  assumption  of  a 
colUsionless  sheath  and  are,  therefore,  valid  only  at  low  pressures  and/or  for 
Ion  species  which  do  not  undergo  significant  resonance  charge  transfer 
collisions.  It  is  anticipated  that  the  incorporation  of  the  effects  of  collisions 
will  be  to  smear  out  the  resulting  ion  bombardment  energy  distribution  and 
displace  the  high  energy  peak  in  the  distribution  toward  lower 
bombardment  energy.  This  result  would  be  consistent  with  previous 
experimental  studies  were  collisions  in  the  sheath  were  an  important 
factor.  5-®  An  extension  of  the  model  to  include  the  effect  of  collisions  is 
currently  in  progress. 

A  further  constraint  on  the  validity  of  these  results  is  imposed  by  the 
assumptions  of  a  constant  electron  temperature  and  charged  particle 
density  in  the  plasma.  As  discussed  in  Ref.  1.  the  energy  modulation  of  the 
majority  of  the  electrons  is  small  for  frequencies  larger  than  approximately 
50  KHz  and  gas  pressures  smaller  than  approximately  1  Torr,  for  the  case  of 
argon  discharges  discussed  in  this  paper.  Even  for  other  gases  (e.g.,  CF4. 
SiH4,  etc.),  this  assumption  is  probably  still  valid  in  the  low  pressure  regime 
required  by  the  collisionless  sheath  assumption.  The  validity  of  the 
assumption  of  a  constant  charged  particle  density  in  the  plasma  volume 
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depends  on  the  loss  rate  of  the  charged  particles  from  the  plasma  and  the 
production  rate  of  the  particles  in  the  plasma.  In  particular,  the  modulation 
in  the  production  rate  (  and,  therefore,  of  the  density)  of  the  charged 
particles  depends  on  the  maintenance  mechanism  of  the  discharge,  as 
discussed  in  Ref.  1 .  The  production  mechanism  of  primary  importance  is, 
however,  a  function  of  the  discharge  parameters  such  as  frequency  of  the 
applied  rf  voltage,  gas  pressure,  gas  mixture,  etc.  The  possible  effect  of  this 
production  rate  modulation  is  not  included  in  the  present  model. 
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SECTION  IV 


AN  EXPERIMENTAL  STUDY  OP  ION  BOMBARDMENT  ENERGY 
DISTRIBUTIONS  IN  RF  DISCHARGES 
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An  Experimental  Study  of  Ion  Bonibardment  Energy  Distributions 

in  RF  Discharges 
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Department  of  Electrical  Engineering 
University  of  Minnesota 
Minneapolis,  MN  55455 

I.  INTRODUCTION 

The  energy  with  which  Ions  In  a  plasma  strike  a  surface  exposed  to 
the  plasma  is  the  subject  of  much  current  Interest.  This  Is  especially  true  for 
radio  frequency,  l.e.,  rf,  driven  plasmas  which  are  used  by  the 
microelectronics  Industry  In  the  fabrication  of  Integrated  circuits.  Two  such 
examples  of  the  use  of  rf  plasmas  are  plasma  etching  and  plasma  deposition. 
In  both  plasma  etching  and  plasma  deposition,  the  energy  with  which  Ions 
strike  the  substrate  surface  can  be  very  critical  In  determining  the 
properties  and  quality  of  the  process.  For  example,  the  Ion  bombardment 
energy  distribution  Is  known  to  effect  the  amount  of  surface  damage  and  the 
degree  of  anisotropy  during  rf  plasma  etching^. 

Measurements  of  the  Ion  bombardment  energy  distribution  of  ions 
Incident  on  an  electrode  in  a  rf  plasma  reactor  have  been  reported 
previously  by  other  researchers.  Cobum  and  coworkers  developed  an  ion 
energy  analysis  system  which  combined  a  90®  deflection  electrostatic 
hemispherical  bandpass  energy  analyzer  with  a  quadrupole  mass  analyzer. 
This  system  was  Initially  described  in  a  study^  of  ion  bombardment  energy 
distributions  In  dc  discharges  used  for  dc  sputter  deposition.  This  system 
was  subsequently  modified  to  measure  the  Ion  bombardment  energy 
distributions  that  occur  In  rf  diode  sputter  deposition^.  In  this  publication, 
these  ion  bombardment  energy  distributions  were  then  utilized  to 
determine  the  maximum  ion  bombardment  energy  for  several  different 
discharge  situations.  However,  only  a  single  set  of  Ion  bombardment  energy 


^D.  L.  Flanun.  V.  M.  Donnelly,  and  D.  E.  Ibbotson,  In  VLSI  Electronics,  edited  by 

N.  G.  Einspnich  and  D.  M.  Brown  (Academic,  Orlando,  1984),  Vol.  8,  Chap.  8,  p.  189. 

2j.  W.  Cobum,  Rev.  ScL  Instrumen.  41,  1219  (1970). 

3j.  W.  Cobum  and  E.  Kiy.  J.  Appl.  Phys.  43, 4965  (1972). 
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distribution  curves  for  Ha*.  H2O+,  and  Eu+  from  a  discharge  in  argon  driven 
with  100  Watts  of  13.56  MHz  rf  power  were  presented  in  this  paper.  Two 
similar  studies'^-  5  by  these  investigators  included  ion  bombardment  energy 
distributions  for:  Ar2+  In  a  50  mTorr,  13.56  MHz  argon  discharge;  ArH+  in  a 
20  mTorr,  13.56  MHz  argon  discharge;  and  ArH+  in  a  100  kHz.  50  mTorr 
argon  discharge. 

Vasile®  utilized  a  system  consisting  of  a  retarding  grid  energy  analyzer 
in  conjunction  with  a  quadrupole  mass  spectrometer  to  measure  the  ion 
bombardment  energy  distributions  of  F2+  and  F+  ions  in  a  rf  discharge  for  a 
mixture  of  helium  and  fluorine.  RF  discharges  in  chlorine  were  studied  by 
Donnelly.  Flamm.  and  Bruce^  using  a  retarding  grid  energy  analyzer  with  a 
quadrupole  mass  spectrometer.  Thompson,  Allen,  Richards,  and  Sawln® 
studied  the  ion  bombsirdment  energy  distributions  in  discharges  of  gases 
containing  various  fluorine  and  chlorine  compounds.  They  also  used  a 
retarding  grid  energy  analyzer,  but  with  a  plate  as  a  detector.  Thus,  their  ion 
energy  distribution  results  were  not  mass  resolved. 

Other  studies  have  measured  the  ion  bombardment  energy  distribution 
incident  on  the  wall  m  contact  with  the  plasma  rather  than  an  electrode.  A 
study  was  performed  by  Vaslle  and  Smollnsky®  using  a  reactant  gas  mixture 
of  argon  and  vlnyltrlmethylsilane.  Measurements  have  also  been  made  of  the 
ion  bombardment  energy  distributions  in  a  rf  driven  Thonemann  ion  source. 
A  paper  by  Erdi®  contains  a  high  resolution  ion  bombardment  energy 
distribution  obtained  by  using  an  electrostatic  deflection  energy  analyzer. 

The  study  described  here  was  undertaken  as  part  of  an  effort  to 
provide  a  database  on  rf  discheirges  for  use  in  the  modeling  of  rf  reactors. 

The  need  for  accurate  data  on  rf  discharges  which  arose  during  the 
development  of  a  model  for  rf  discharges  by  Metze,  Ernie,  and  Oskami* 


^K.  Kohler,  J.  W.  Cobum.  D.  E.  Home,  E.  Kay.  andJ.  H.  Keller,  J.  Appl.  Phys.  57,  59  (1985). 

®K.  KOhler,  D.  E.  Home,  and  J.  W.  Cobum,  J.  Appl.  Phys.  58,  3350  (1985). 

®M.  J.  Vaslle.  J.  Appl.  Phys.  51.  2503  (1980). 

^V.  M.  Donnelly  and  D.  L.  Flamm.  J.  Appl.  Phys.  58.  2135  (1985). 

®B.  E.  Thompson.  K.  D.  Allen.  A.  D.  Richards,  and  H.  H.  Sawln.  J.  Appl.  Ph3rs.  59.  1890  (1986). 
®M.  J.  Vaslle  and  G.  Smollnsky.  Int.  J.  Mass  Spectrom.  Ion  Phys.  12,  133  (1973). 
lOj.  ErO.  Nuclear  Instim.  3,  303  (1958). 

*  -A.  Metze.  D.  W.  Ernie,  and  H.  J.  Oskam,  J.  Appl.  Ph)^.  60,  3081  (1986). 


35 


demonstrated  the  necessity  for  this  experimental  work.  This  model  has 
been  used  to  predict  the  ion  bombardment  energy  distributions  at  the 
electrodes  of  a  capacitively  coupled  rf  discharge^^.  Since  only  limited 
experimental  data  for  the  ion  bombardment  energy  distribution  in  rf 
discharges  was  available,  these  experiments  were  undertaken  to  increase 
the  range  of  pressures  and  frequencies  over  which  the  ion  bombardment 
energy  distributions  have  been  measured  so  that  a  more  realistic  comparison 
can  be  made  Avith  the  model  predictions.  Since  it  is  desirable  to  first 
understand  the  physics  of  the  rf  plasma,  noble  gases  were  used  in  these 
experiments  in  order  to  simplify  complications  due  to  chemical  reactions  in 
the  plasma.  It  is  hoped  that  the  refinement  of  present  models  of  a  rf  plasma 
using  the  results  of  these  experiments  will  ultimately  lead  to  a  better 
understanding  of  more  complicated  plasma  systems,  such  as  those  currently 
used  in  the  plasma  processing  environment. 


n.  EXPERIBAENTAL  APPARATUS 


A.  General  Configuration 

The  vacuum  system  used  in  these  experiments  is  shown  schematically 
in  Fig.  1 .  This  vacuum  system  consisted  of  two  separate  chambers.  The  first 
chamber  is  the  reactor  chamber  in  which  the  desired  rf  discharge  was 
maintained.  The  second  chamber  is  the  analysis  chamber  which  housed  a 
retarding  grid  energy  analyzer  and  a  quadrupole  mass  spectrometer.  The 
two  chambers  were  separated  by  a  wall  containing  a  100  nm  diameter 
sampling  orifice  which  allowed  for  sampling  of  the  particles  that  were 
Incident  on  this  wall  of  the  reactor  chamber  followed  by  subsequent  analysis 
of  these  sampled  particles  in  the  analysis  chamber. 


12a.  Metze.  D.  W.  Ernie,  and  H.  J.  Oskam,  37*  Gaseous  Electronics  Conference,  Bull.  Amer. 
Phys.  Soc.  30.  145  (1985). 
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FIG^^^^lock^dlagran^^h^j^cuun^^sten^ 


A  high  vacuum  valve  Interconnected  the  reactor  chamber  and  the 
analysis  chamber  and  wss  used  for  evacuating  the  reactor  by  the  analysis 
chamber’s  vacuum  pumps  when  the  system  wss  not  in  use.  During  data 
acquisition,  this  valve  was  closed  £ind  the  plasma  reactor  chamber  was 
maintained  at  a  constant  pressure  by  an  electronic  feedback  control  system 
which  used  a  capacitance  manometer  as  input  and  controled  an  electrically 
controlled  Inlet  valve  which  was  connected  to  the  gas  source.  A  two  stage 
rotary  vane  pump  was  also  connected  to  the  reactor  chamber  through  a 
second  high  vacuum  valve  and  a  butterfly  valve.  The  system  was  normally 
operated  In  a  quasi-static  mode  with  these  two  valves  closed.  In  the  quasi- 
static  mode,  the  only  loss  of  gas  from  the  reactor  was  through  the  sampling 
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orifice.  In  order  to  compensate  for  this  loss,  a  equal  amount  of  reactant  gas 
was  introduced  Into  the  plasma  reactor  through  the  electrically  controlled 
inlet  valve.  By  this  means,  the  pressure  in  the  reactor  was  held  at  a  constant 
value.  In  this  mode  of  operation,  the  mechanical  pump  was  simply  used  to 
rough  down  the  reactor  chamber  before  interconnecting  it  with  the  analysis 
chamber  for  further  evacuation.  In  addition,  however,  the  butterfly  valve 
allowed  the  system  to  also  be  used  as  a  flowing  system  where  the  pressure 
was  again  maintained  at  a  constant  value,  but  the  gas  loss  was  than 
determined  by  both  the  loss  due  to  the  orifice  and  the  loss  through  the 
butterfly  valve.  In  this  mode,  the  mechanical  pump  had  an  active  role  and 
the  flow  rate  was  determined  by  the  butterfly  valve  setting.  In  either  mode, 
the  pressure  control  sjrstem  made  it  possible  to  maintain  any  desired 
constant  pressure  from  10  mTorr  to  10  Torr. 

The  analysis  chamber  housed  the  energy  analyzer,  a  quadrupole  mass 
spectrometer,  and  an  ion  multiplier  detector.  This  chamber  was  evacuated 
by  a  turbomolecular  pump  vdiich  was  backed  by  a  two-stage  rotary  vane 
pump.  When  the  reactor  chamber  was  evacuated,  the  anal3rsis  chamber 
achieved  a  pressure  of  10‘®  Torr.  At  all  pressures  of  the  reactor  chamber 
used  in  the  experiments,  the  pressure  of  the  analysis  chamber  was  always 
below  10-5  Torr.  It  was  necessary  to  maintain  this  low  pressure  in  the 
analysis  chamber  so  that  the  sampled  particles  from  the  reactor  undergo  no 
additional  collisions  in  the  analysis  chamber.  The  pressure  in  this  chamber 
was  monitored  by  an  ionization  gauge. 

The  energy  analyzers  used  were  of  the  retarding  grid  type.  A  retarding 
grid  energy  analyzer  utilizes  a  metal  mesh  or  grid  vdiich  is  attached  to  a 
variable  voltage  source.  A  beam  of  charged  particles  with  a  certain  kinetic 
energjr  directed  towsirds  the  grid  can  only  reach  and  pass  through  the  grid  if 
the  kinetic  energy  of  the  particles  is  greater  than  the  product  of  the  charge 
on  the  particle  times  the  electric  potential  of  the  grid  measured  with 
respect  to  the  electric  potential  at  the  location  vdiere  the  particles 
originated.  A  properly  designed  retarding  grid  energy  Einalyzer  allows  all  the 
particles  with  a  kinetic  energy  to  charge  ratio  greater  than  the  retarding 
voltage  to  pass  through  the  grid  and  be  detected  by  the  detection  system. 
Thus,  it  should  be  noted  that  a  retarding  grid  energy  analyzer  yields  an 
integral  of  the  ion  energy  probability  density  function  and  not  the  function 
itself.  However,  the  ion  energy  probability  density  function  may  be  found  by 
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taking  the  derivative  of  the  measured  result  with  respect  to  the  ion  energy. 
Actually,  for  many  uses,  one  is  more  interested  in  the  niamber  of  ions  having 
an  impact  energy  greater  (or  less)  than  a  certain  energy.  This  is  the 
parameter  that  was  directly  measured  by  a  retarding  grid  energy  analyzer. 

The  quadrupole  mass  analysis  ^tem  consisted  of  a  quadrupole  mass 
filter  and  an  electrostatic  lense  array  to  help  focus  the  incoming  ion  beam 
into  the  quadrupole.  There  was  also  an  electron  impact  ionization  assembly 
within  the  lense  array.  This  allowed  mass  analysis  of  the  neutr2d  content  of 
the  beam  extracted  through  the  orifice  from  the  plasma.  The  ion  multiplier 
detector  was  run  in  a  pulse  coimting  mode  in  order  to  achieve  high 
sensitivity.  The  resulting  mass  analysis  system  had  enough  sensitivity  to 
allow  for  a  mass  analysis  of  residual  background  gases  when  the  entire 
system  was  evacuated  to  10'®  Torr.  Besides  the  vacuum  system,  there  were 
support  electronics  for  maintaining  and  biasing  the  plasma  and  for  the 
quadrupole  mass  spectrometer  and  its  ion  multiplier  detector  system. 
Further  support  equipment  was  also  provided  for  the  energy  anal3rsis  system. 
A  computer  controlled  data  acquisition  system  was  used  to  coordinate  the 
operation  of  the  quadrupole  and  energy  analysis  systems. 

B.  Initial  Design  of  Reactor  and  Energy  Analyzer 

Two  different  reactor  chamber  designs  were  used.  The  first  reactor 
was  of  the  cylindrical  discharge  tube  type  and  is  diagrammed  in  Fig.  2.  The 
discharge  tube  was  constructed  of  2.5  cm  diameter  I^ex  tubing.  The 
sampling  orifice  was  located  in  a  narrow  metal  band  which  replaced  a  small 
section  of  the  Pyrex  tubing.  A  100  pm  diameter  orifice  was  machined  in  the 
band  at  a  location  where  the  band  was  20  pm  thick.  This  combination  of 
metal  thickness  and  orifice  size  was  chosen  in  order  to  minimize  collisions 
with  the  sides  of  the  orifice  and  to  ensure  that  the  flow  through  the  orifice 
was  in  a  molecular  regime.  The  discharge  tube  contained  three  additional 
electrodes,  any  given  pair  of  which  could  be  chosen  as  the  anode  and 
cathode.  These  electrodes  were  initially  short  sections  of  hollow  stainless 
steel  tubing.  However,  other  electrode  shapes  for  the  cathode  were  tried  in 
order  to  alleviate  problems  with  noise  in  the  discharge,  as  discused  below. 
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To  Gas  Inlet  Valve 


In  the  experiments  which  utilized  the  cylindrical  reactor  shown  in 
ng.  2,  the  reactor  was  biased  as  shown  in  Fig.  3.  The  dc  current  in  the 
discharge  tube  was  set  by  supply  Vq.  The  dc  potential  of  the  anode  with 
respect  to  ground  was  determined  by  supply  Vb.  As  a  result.  Vb  could 
effectively  be  used  to  appty  a  dc  bias  potential  to  the  sampling  electrode 
with  respect  to  the  potential  of  the  plasma  at  the  location  of  the  sampling 
electrode.  This  was  done  by  first  allowing  the  sampling  electrode  to  float 
and  measuring  the  voltage  on  the  sampling  electrode  with  respect  to 
ground,  yielding  the  sum  of  the  electric  potential  of  the  plasma  at  the 
location  of  the  sampling  electrode  and  the  floating  potential  of  the  wall  with 
respect  to  the  plasma.  This  voltage  could  be  adjusted  by  adjusting  Vb.  Then, 
when  the  sampling  electrode  was  subsequently  grounded,  a  voltage  equ£d  to 
the  floating  potential  minus  the  previously  measured  voltage  was  effectively 
applied  between  the  plasma  and  the  sampling  electrode. 
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Supply  Vrf  provided  an  additional  small  rf  potential  to  the  discharge 
region  without  significantly  disturbing  the  dc  conditions  of  the  tube.  It 
should  be  noted  that  even  though  the  anode  and  cathode  of  the  discharge 
tube  were  at  different  dc  potentials,  they  (and  the  entire  plasma  In  the 
discharge  tube)  were  at  the  same  rf  potential.  This  biasing  arrangement 
allowed  for  a  known  rf  potential  superimposed  with  an  adjustable  dc 
potential  to  be  applied  between  the  plasma  and  the  sampling  electrode.  This 
biasing  arrangement  corresponded  to  the  conditions  across  a  single  sheath 
used  In  developing  the  Ion  bombardment  energy  distribution  model 
previously  mentioned,  l.e..  Refs.  11  and  12.  It  was  hoped  that  this  setup 
could  yield  Information  on  the  effect  of  a  time  varying  sheath  potential  on 
the  energy  distribution  of  Ions  traversing  the  sheath. 


FIG.  3.  Biasing  scheme  used  with  the  cylindrical  reactor. 


The  initial  experimental  setup  utilized  the  cylindrical  discharge  tube 
described  above  and  a  three  grid  planar  retarding  grid  energy  analyzer.  The 
energy  analyzer  was  located  at  the  entrance  to  the  focusing  lenses  of  the 
quadrupole  mass  spectrometer.  This  setup  failed  to  3rleld  acceptable  results 
for  several  reasons.  First,  the  electrostatic  lenses  of  the  quadrupole  are 
themselves  energy  selective.  A  small  variation  in  the  bias  potential  of  the 
electrostatic  lenses  produced  large  changes  in  the  resulting  ion  energy 
distribution  measurement.  An  attempt  to  correct  this  problem  was  made  by 
removing  the  electrostatic  lenses  from  the  quadrupole  and  mounting  the 
energy  analyzer  directly  on  the  quadrupole  Itself.  However,  this  failed  to 
yield  the  desired  Improvement  because  the  quadrupole  itself  is  also  energy 
selective.  However,  an  additional  and  fer  worse  problem  was  the  criticality  of 
the  alignment  of  the  quadrupole  and  the  orifice,  with  small  changes  in  the 
alignment  of  the  quadrupole  relative  to  the  orifice  resulting  in  large 
deviations  in  the  measured  ion  energy  distribution.  It  was  realized  that  if  the 
incoming  ion  beam  is  not  perfectly  perpendicular  to  the  energy  analyzer,  the 
varying  potential  on  the  second  grid  would  cause  the  beam  to  be  deflected 
by  varjdng  amounts  depending  upon  the  potential  applied  to  this  grid  and 
the  velocity  of  the  individual  ions.  This,  in  turn,  caused  the  beam  to  be 
misaligned  by  varying  amoimts  when  entering  the  quadrupole  and,  hence, 
the  amount  of  the  beam  actually  detected  had  an  additional  functional 
dependence  on  of  the  potential  applied  to  the  retarding  grid.  In  fact,  it  was 
possible  to  align  the  system  so  that  the  amount  of  signal  detected  increased 
over  a  region  where  the  retarding  potential  for  positive  ions  was  also 
increasing.  This  is  clearly  an  invalid  result,  since  an  increase  in  the 
retarding  potential  can  only  result  in  a  decrease  In  the  number  of  positive 
Ions  passing  through  the  energy  analyzer.  This  increase  In  signal  was  not  due 
to  an  increase  in  the  number  of  positive  ions  transmitted  through  the  energy 
analyzer,  but  due  to  the  Ion  beam  being  swept  Into  alignment  with  the 
quadrupole  by  the  energy  analyzer  so  that  the  transmission  losses  of  the 
quadrupole  were  reduced.  Because  of  these  problems  it  was  decided  that, 
although  the  quadrupole  could  be  used  to  determine  which  ionic  species 
were  present,  it  could  not  be  used  in  conjunction  with  a  retarding  grid 
energy  analyzer  to  yield  unambiguous  ion  bombardment  energy  distributions. 
Subsequent  experiments  were  performed  using  a  plate  (or  mesh)  behind  the 
energy  analyzer  as  the  detector.  While  this  scheme  did  not  allow  for  the 
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measurement  of  mass  analyzed  ion  bombardment  energy  distributions,  it  did 
yield  reliable  results.  Furthermore,  by  choosing  the  correct  discharge 
conditions,  the  sampled  ion  beam  could  be  made  to  consist  of 
predominantly  one  ion  species,  as  confirmed  by  mass  spectrometer  analysis. 

This  resulting  setup  still  failed  to  yield  the  desired  results  due  to 
oscillations  in  the  dc  discharge.  It  was  found  that  when  there  was  no  applied 
rf  potential  and  the  sampling  electrode  was  allowed  to  float,  the  voltage  of 
the  sampling  electrode  tended  to  oscillate.  This  oscillation  generally  had  an 
amplitude  of  about  ten  volts  and  was  strongly  affected  by  the  current  through 
the  tube  and  the  gas  pressure  within  the  tube.  Since  previous  studies ^3  had 
Indicated  that  the  noise  in  dc  discharges  may  be  a  strong  function  of  anode 
size  and  decreases  for  Increasing  anode  area,  attempts  were  made  to 
eliminate  these  oscillations  by  increasing  the  area  of  the  anode  and  by  using 
a  planar  anode  surface.  While  these  new  design  features  did  decrease  the 
oscillations  slightly,  this  decrease  was  not  sufficient  to  allow  the 
measurements  to  be  performed  at  the  required  pressures.  Because  of  this 
oscillation  problem,  the  cylindrical  discharge  tube  design  was  abandoned 
and  the  planar  reactor  described  in  the  next  section  was  constructed. 

C.  Final  Design  of  Reactor  and  Energy  Analyzer 

The  second  reactor  chamber  was  a  parallel  plate  reactor  as  shown  in 
Fig.  4.  The  reactor  chamber  consisted  of  a  10  cm  I^ex  bell  Jar  mounted  on 
a  15  cm  adapter  plate  in  the  center  of  which  was  a  100  ^m  diameter  orifice 
vdiich  connected  the  reactor  chamber  to  the  analysis  chamber.  The  adapter 
flange  served  as  one  electrode  for  the  discharge.  The  second  electrode  was  a 
9  cm  metal  plate  which  ccould  be  placed  a  variable  distance  from  the  lower 
electrode.  This  reactor  realistically  represented  the  geometry  for  which  the 
equivalent  circuit  model  of  an  rf  discharge  given  in  Ref.  1 1  was  derived. 


^3a.  Garacadden.  in  Caseous  Electronics,  edited  by  M.  N.  Hirsh  and  H.  J.  Oskam  (Academic. 
New  York.  1978).  Vol.  1.  Chap.  2.2.  p.  65. 
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FIG.  4.  Design  of  the  parallel  plate  reactor  and  the  hemispherical  energy 
analyzer _ 


In  the  series  of  experiments  described  in  section  III  which  utilized 
the  planar  reactor  shown  in  Fig.  4.  the  planar  reactor  was  driven  by  a 
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variable  frequency  rf  power  supply.  The  method  of  applying  the  rf  power  to 
the  parallel  plate  reactor  is  shown  In  Fig.  5.  It  consisted  of  the  power  supply 
and  a  matching  network  which  was  connected  through  a  blocking  capacitor, 
Cb,  to  the  driven  electrode.  Different  matching  networks  were  used, 
depending  upon  the  frequency  of  operation  desired.  One  successful  method 
was  to  use  a  transformer  woimd  on  a  ferrite  core  and  designed  for  the 
desired  frequency  of  operation.  The  use  of  the  blocking  capacitor  was  most 
useful  at  higher  frequencies,  since  it  then  allows  the  time-averaged  (dc) 
current  through  the  reactor  to  be  measured  by  grounding  the  blocking 
inductor  shown  in  Fig.  5  through  a  dc  ammeter.  In  the  same  manner,  the 
average  dc  potential  at  the  driven  electrode  can  be  determined  by 
connecting  the  blocking  inductor  to  a  dc  voltmeter.  A  variation  of  this 
method  was  used  at  lower  frequencies,  where  the  blocking  capacitor  Cb,  the 
blocking  inductor,  and  the  ammeter  and  voltmeter  were  omitted. 


Matching  Blocking 

Network  Inductor 


FIG^^S^Biasln^iet^rit^use^^d^h^par^lelplatereartor. 


In  order  to  increase  the  resolution  of  the  energy  analyzer  and  to 
alleviate  the  alignment  problems  previously  mentioned,  the  energy  analyzer 
also  underwent  a  major  design  revision.  This  final  version  of  energy  anal3^er 
utilized  three  hemispherical  grids  with  a  fourth  hemispherical  plate  (or 
grid)  being  used  as  a  detector.  Since  the  grids  were  hemispherical  and  the 
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orifice  was  located  at  the  center  of  curvature  of  the  grids,  the  retarding 
electric  field  Is  parallel  to  the  Ion’s  velocity.  This  will  help  prevent  any  loss 
of  resolution  by  limiting  the  deflection  of  the  Ions  from  their  original  path. 
Since  the  detector  was  also  hemispherical,  the  detection  of  the  ions  no 
longer  depends  critically  on  the  alignment  of  the  orifice  and  the  detector. 
Initially  the  detector  was  a  grid.  This  grid  was  subsequently  replaced  by  a 
plate  in  order  improve  the  signal  sensitivity  by  collecting  all  of  the  ion  beam 
rather  than  just  the  fraction  that  was  collected  by  the  partially  transparent 
grid.  As  expected,  no  major  changes  In  the  shape  of  the  measured  ion 
energy  distribution  function  were  found  vdien  changing  from  the  grid 
detector  to  the  plate  detector.  However,  as  expected,  this  resulted  in  a  large 
increase  In  the  detected  Ion  current. 

Considering  the  grids  shown  In  Fig.  4  to  be  numbered  according  to 
their  order  from  the  orifice,  the  first  grid,  the  orifice,  and  the  vacuum 
housing  were  all  at  groimd  potential  so  that  the  ion  trajectories  were 
undisturbed  until  the  ions  passed  through  the  first  grid.  The  second  grid's 
voltage  was  controlled  by  the  data  acquisition  system  and  was  varied  to 
exclude  all  ions  with  a  kinetic  energy  to  charge  ratio  less  than  a  desired 
value.  The  third  grid  was  placed  at  a  fixed  potential  of  -5  Volts  In  order  to 
block  electrons  which  had  also  been  sampled  through  the  orifice  and  to 
repel  any  secondary  electrons  emitted  from  the  detector  back  onto  the 
detector.  The  signal  from  the  detector  plate  was  interfaced  to  the  data 
acquisition  system. 

The  operation  of  the  planar  plasma  reactor  and  the  hemispherical 
retarding  grid  energy  analyzer  was  verified  e^erimentally.  This  was  done  by 
measuring  the  energy  distribution  of  positive  ions  from  a  dc  discharge 
where  the  upper  electrode  was  used  as  the  anode  and  the  lower  electrode 
containing  the  sampling  orifice  was  used  as  the  cathode.  The  results  found 
in  these  measurements  were  consistent  with  the  results  reported  previously 
by  Davis  and  Vandersllcei^.  Subsequently,  experiments  were  performed 
using  rf  voltages  across  the  reactor.  The  results  of  these  experiments  will  be 
described  in  the  following  sections.  The  resialts  and  behavior  of  the 
experimental  setup  under  edl  these  conditions  indicate  that  the  system 
performed  as  desired. 


^^W.  D.  Davis  and  T.  A.  Vanderslice,  Phys.  Rev.  131. 219  (1963). 
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m.  EXPERIMENTAL  METHOD  AND  RESULTS 


Using  the  final  version  of  the  experimental  setup  described  in  Section 
II.  various  experiments  were  performed  in  which  a  single  experimental 
parameter  was  varied  while  the  other  parameters  were  held  constant.  The 
independent  parameters  considered  in  these  experiments  were:  the 
pressure  of  the  reactant  gas  in  the  reactor,  the  tsqje  of  reactant  gas  (or 
reactant  gas  mixture),  the  area  ratio  of  the  electrode  surfaces,  the  electrode 
spacing,  the  magnitude  of  the  voltage  across  the  reactor,  and  the  frequency 
of  the  voltage  across  the  reactor.  Results  on  the  study  of  the  ion 
bombardment  energy  distribution  using  this  experimental  setup  are 
presented  below,  grouped  according  to  the  discharge  parameter  varied. 


A.  Effect  of  Varying  Pressure 

Data  was  obtained  with  a  13.6  MHz,  200  V  peak-to^eak  voltage 
across  the  reactor  for  different  values  of  gas  pressure  with  neon  as  the 
reactant  gas.  Figure  6a  shows  the  ion  current  collected  as  a  fimctlon  of  the 
retarding  grid  voltage.  The  ion  bombardment  energy  distribution  is  shown  in 
Fig.  6b  and  is  the  derivative  of  the  data  shown  in  Fig.  6a.  As  shown  in  Fig.  6b, 
for  a  pressure  of  25  and  50  mTorr,  the  energy  distribution  of  the  neon  ions 
is  basically  a  single  narrow  peak  at  an  energy  of  50  eV.  The  distribution 
drops  off  rapidly  from  its  peak  value  for  energies  greater  than  50  eV. 
However,  the  distribution  drops  off  more  slovdy  and  in  a  manner  that 
appears  exponential  for  energies  less  than  50  eV.  However,  at  a  pressure  of 
100  mTorr,  the  low  energy  side  of  the  peak  drops  off  much  slower,  so  that 
the  distribution  decreases  almost  linearly  with  decreasing  energy.  The  high 
energy  side  of  the  peak  remains  relatively  unchanged.  At  200  mTorr,  the 
distribution  takes  on  an  almost  symmetric  shape  with  the  peak  in  the 
distribution  now  shifted  to  the  center  of  the  distribution.  At  400  and  800 
mTorr,  the  peak  shifts  towards  even  lower  energies.  In  addition,  the 
maximum  detected  ion  bombardment  energy  also  decreases  as  the  pressure 
increases  for  these  curves. 
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FIG.  6a.  Hie  ion  current  detected  as  a  function  of  Ion  bombardment  energy. 
The  data  was  obtained  with  a  13.6  MHz  200  V  peak-to^eak  sine  wave 
applied  to  the  upper  electrode.  The  pressure  of  neon  In  the  reactor  was:  (a) 
25  mTorr;  (b)  50  mTorr;  (c):  100  mTorr;  (d)  200  mTorr;  (e)  400  mTorr;  (f) 
800  mTorr. 
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FIG.  6b.  The  ion  bombardment  energy  distribution  calculated  from  the  data 
shown  in  Fig.  6a.  The  data  was  obtained  using  a  13.6  MHz  200  V  peak-to- 
peak  sine  wave  applied  to  the  upper  electrode.  The  pressure  of  neon  in  the 
reactor  was:  (a)  25  mTorr;  (b)  50  mTorr;  (c):  100  mTorr;  (d)  200  mTorr; 

(e)  400  mTorr;  (f)  800  mTorr. 

The  changes  in  the  ion  bombardment  energy  distribution  as  a  result  of 
the  changes  in  the  pressure  of  neon  in  the  reactor  chamber  are  consistent 
with  vdiat  is  expected.  At  pressures  below  50  mTorr,  the  ions  that  traverse 
the  sheath  do  so  without  a  significant  number  of  these  ions  experiencing  a 
collision.  This  is  expected  since  the  mean  free  path  of  a  neon  ion  in  neon  is 
approximately  1  cm  at  a  pressure  of  25  mTorr,  while  the  sheath  has  a 
dimension  on  the  order  of  1  mm.  As  the  pressure  increases,  the  mean  free 
path  of  the  ion  decreases  and  at  200  mTorr  the  mean  free  path  is 
approximately  1  mm.  Therefore,  at  200  mTorr,  it  expected  that  only  about 
half  of  the  ions  traverse  the  sheath  without  experiencing  a  collision.  Since  in 
neon,  the  predominant  collision  mechanism  is  symmetric  resonant  charge 
transfer,  ^  5  the  ions  that  undergo  collision  will  transfer  almost  all  of  their 

Chapman,  Gtotu  Dlsdiarge  Processes  (Wiley.  New  York.  1980),  p.39. 
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kinetic  energy  to  the  neutral  neon  atoms.  Thus  the  energy  with  which  an  ion 
bombards  the  surface  will  be  the  energy  that  the  ion  gains  from  the  electric 
field  in  the  sheath  after  its  last  collision  before  reaching  the  surface.  Since  at 
200  mTorr,  most  ions  undergo  a  single  collision,  the  peak  in  the  energy 
distribution  is  shifted  to  half  of  the  energy  of  the  peak  at  25  mTorr.  At  400 
mTorr,  the  mean  free  path  for  an  ions  is  0.5  mm  and  most  ions  undergo  a 
collision  a  fourth  of  the  sheath  distance  from  the  surface  and  the  energy  of 
the  peak  in  the  distribution  is  shifted  to  a  fourth  of  the  25  mTorr  value. 
Similarly,  at  800  mTorr,  the  mean  free  path  is  0.25  mm  and  most  ions 
undergo  a  collision  an  eighth  of  the  sheath  distance  and  the  energy  of  the 
peak  Is  roughly  an  eighth  of  the  25  mTorr  value.  This  analysis  is  only 
approximate  since  it  is  known  that  the  variation  of  the  potential  in  the 
sheath  is  not  linear  and  thus  the  collision  distances  used  above  are  only 
rough  estimates.  In  addition,  the  sheath  dimension  is  actually  a  function  of 
the  pressure.  However,  this  simple  analysis  does  show  that  the  measured 
effect  of  pressure  on  the  ion  bombardment  energy  distribution  is  consistent 
with  known  data  about  the  mean  free  path  of  neon  ions  in  neon. 

One  further  interesting  observation  is  that  the  total  ion  current 
collected  does  not  depend  upon  the  pressure  from  25  mTorr  to  200  mTorr 
for  these  measurements  at  13.6  MHz  with  the  voltage  across  the  reactor 
held  constant  at  200  V  peak-peak.  This  would  imply  that  the  density  of  ions 
and  electrons  is  independent  of  the  gas  pressure  in  this  range  and  is  only 
dependent  on  the  reactor  voltage. 

B.  Effect  of  Varying  Frequency 

Observations  on  the  effects  of  varying  frequency  were  performed  in 
neon  at  pressures  from  50  to  a  500  mTorr.  While  this  work  was  only 
qualitative  in  nature,  one  observation  of  relevance  Is  that  as  the  frequency  is 
lowered,  the  voltage  required  across  the  reactor  to  maintain  a  discharge 
increases  dramatically.  For  example,  below  a  few  MHz,  the  peak-to-peak 
voltage  required  to  sustain  the  discharge  Increases  to  several  hundred  Volts, 
while  at  13.6  MHz  a  discharge  may  be  maintained  with  a  peak-to-peak 
voltage  of  less  than  50  Volts.  These  observations  show  that  the  sustaining 
mechanisms  of  the  plasma  are  influenced  by  the  frequency  at  which  the 
plasma  is  driven.  At  lower  frequencies  (below  a  few  hundred  kHz),  when  an 
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electrode  is  at  a  negative  potential  with  respect  to  the  other  electrode,  it 
will  act  in  a  manner  which  is  similar  to  a  cathode  In  a  dc  discharge. 
Therefore,  the  production  of  ion-electron  pairs  near  the  electrode  will  be 
the  result  of  electron  Impact  ionization  due  to  secondauy  electrons  which 
sure  emitted  by  ion  bombardment  of  the  electrode  and  which  are  accelerated 
back  into  the  plasma.  It  should  be  noted  that  this  mechanism  is  also  present 
at  higher  frequencies  (above  a  few  MHz).  In  addition,  at  lower  frequencies, 
the  electron  temperature  and  possibly  the  electron  density  will  be  time 
varying  and  modulated  by  the  driving  rf  voltage.  At  higher  frequencies,  the 
situation  is  different.  At  these  frequencies,  the  electron  temperature  in  the 
plasma  will  be  time  independent  and  the  high-energy  electrons  in  the  tall  of 
the  electron  energy  distribution  will  contribute  significantly  to  the 
production  of  ion-electron  palrs^®.  This  contribution  is  less  important  at  the 
lower  frequencies  since  the  electrons  do  not  reach  equilibrium. 

C.  Effect  of  Vaiying  Voltage 

The  data  in  Fig.  7  was  obtained  from  a  rf  discharge  in  neon  at  a 
frequency  of  13.6  MHz  and  at  a  pressure  of  50  mTorr.  Figure  7a  shows  the 
measured  ion  current  and  Fig.  7h  shov^  the  resulting  ion  bombardment 
energy  distribution.  The  voltage  across  the  reactor  was  a  sine  wave  with  a 
peak-to-peak  amplitude  which  was  varied  from  40  V  to  500  V.  Several 
trends  are  indicated  by  the  data  shown  In  Fig.  7b.  First,  the  width  of  the 
peak  in  the  Ion  bombardment  energy  distribution  appears  unefiected  by 
changes  In  the  magnitude  of  the  voltage  across  the  reactor.  This  peak  width 
is  approximately  15  eV.  However,  the  energy  of  the  peak  increases  as  the 
voltage  across  the  reactor  Increases.  Finally,  the  positive  ion  current 
detected  also  increases  as  the  voltage  across  the  reactor  increases. 


16m.  GIU.  Vacuvun  34.  357  (1984). 
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Ion  Bombardment  Energyr  (eV) 

FIG.  7a.  The  ion  current  detected  as  a  function  of  ion  bombardment  energy. 
The  pressure  of  neon  in  the  reactor  was  50  mTorr.  A  13.6  MHz  sine  wave 
was  applied  to  the  upper  electrode  with  a  peak-to-peak  voltage  of:  (a)  500 
V;  (b)  400  V;  (c)  300  V;  (d)  200  V;  (e)  160  V;  (f)  120  V;  (g)  80  V;  (h)  40  V. 
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FIG.  7b.  The  ion  bombardment  energy  distribution  as  calculated  from  the 
data  presented  in  Fig.  7a.  The  pressure  of  neon  in  the  reactor  was  50  mTorr. 
A  13.6  MHz  sine  wave  was  applied  to  the  upper  electrode  with  a  peak-to- 
peak  voltage  of:  (a)  500  V;  (b)  400  V;  (c)  300  V;  (d)  200  V;  (e)  160  V;  (f) 

120  V;  (g)  80  V:  (h)  40  V. 

In  Fig.  8,  the  peak  in  the  ion  bombardment  energy  for  the 
distributions  given  in  Fig.  7b  are  graphed  as  a  function  of  the  peak-to-peak 
voltage  across  the  reactor.  This  graph  shows  that  the  peak  Ion  bombardment 
energy  Epeak  Is  a  linear  function  of  the  peak-to-peak  voltage  Vp.p  across  the 
reactor.  Using  the  model  presented  In  Refs.  11  and  12,  one  can  predict  the 
dependence  of  the  peak  in  the  ion  bombardment  energy  on  the  pe8ik-to- 
peak  voltage  across  the  reactor.  For  a  planar  rf  reactor  with  approximately 
equal  area  electrodes,  the  voltage  across  the  sheath  at  an  electrode  can  be 
approximated  by  the  floating  potential  during  half  the  rf  cycle.  During  the 
other  half  of  the  rf  cycle,  the  voltage  may  be  approximated  by  the  floating 
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potential  plus  the  positive  half  of  the  applied  sine  wave.  Since,  at 
frequencies  in  the  MHz  range,  the  ions  take  several  rf  cycles  to  traverse  the 
sheath,  the  ions  which  do  not  undergo  collisions  bombard  the  surface  with 
an  energy  corresponding  to  the  average  voltage  across  the  sheath  times  the 
charge  on  the  ion.  Collisions  would  have  the  effect  of  broadening  the  ion 
bombardment  energy  distribution  around  this  predicted  peak  value.  Under 
the  simplifying  assumptions  given,  the  time  averaged  voltage  across  the 
sheath  is  the  floating  potential  Vf  plus  the  peak-to-peak  voltage  divided  by 
2k.  Thus,  the  peak  ion  bombardment  energy  should  be  given  by 
Epeak/q  =  Vf+  0.159  Vp.p,  where  q  is  the  the  charge  of  the  bombarding 
positive  ion.  However,  as  indicated  by  a  least  squares  fit,  the  data  presented 
in  Fig.  8  is  linear  with  a  slope  of  0.21  rather  than  0.159.  One  possible 
explanation  for  this  discrepancy  is  the  assumption  that  the  ions  have  an 
average  energy  equal  to  the  average  voltage  across  the  sheath  times  the 
charge  on  the  ion.  This  assumption  is  partially  based  on  the  idea  that  the 
ions  enter  the  sheath  with  an  equal  probability  for  all  possible  phase  angles 
of  the  applied  potential.  This  is  not  completely  accurate,  since  the  sheath 
will  change  its  thickness  for  changes  in  the  electric  potential  across  the 
sheath.  Therefore,  it  is  more  probable  that  ions  wUl  enter  the  sheath  during 
that  portion  of  the  rf  cycle  where  the  voltage  across  the  sheath  and,  hence, 
the  sheath  thickness  is  largest.  This  would  lead  to  a  larger  value  for  the  peak 
ion  bombardment  energy,  as  is  observed.  The  intercept  of  the  data 
presented  in  Fig.  8  predicts  a  floating  potential  of  9.3  volts  which  is 
consistent  with  the  predictions  of  the  model  presented  in  Ref.  1 1 . 


D.  Effects  of  Chaxijflxig  the  Electrode  Spacing  and  Area  Ratio 


The  measurements  show  very  little  effect  of  the  electrode  spacing  on 
the  discharge  characteristics  or  the  resulting  ion  bomb8u:dment  energy 
distribution  at  13.6  MHz  and  spacings  between  6  cm  and  15  cm.  Results  at 
13.6  MHz  also  show  little  effect  of  the  area  ratio  on  the  shape  of  the  ion 
bombardment  energy  distribution.  However,  the  time-averaged  dc  potential 
of  the  driven  electrode  is  more  negative  when  the  area  of  the  driven 
electrode  is  reduced.  In  addition,  the  magnitude  of  the  time  averaged  dc 
current  through  the  reactor  when  the  upper  electrode  is  dc  groimded 


through  the  blocking  inductor  also  Increases  as  the  area  of  the  driven 
electrode  is  reduced. 

IV.  IMPUCATIONS  OP  EXPERIBflENTAL  RESULTS  FOR  MODELING  RP 

DISCHARGES 

From  the  data  that  has  been  collected,  one  observation  is  apparent. 

The  effect  of  collisions  in  the  sheath  region  must  be  considered  in  modeling 
the  ion  bombardment  energy  distribution  of  the  positive  ions.  At  pressures 
of  100  mTorr  and  above  in  neon,  the  effect  of  collisions  on  the  measured  ion 
bombardment  energy  distribution  is  quite  dramatic.  In  a  weakly  ionized 
plasma,  the  predominant  collisions  will  be  between  ions  and  neutrals.  Of  the 
different  ion-neutral  collision  mechanisms,  the  most  probable  is 
symmetrical  resonant  charge  transfer  for  helium  and  neon  ions  moving  in 
their  parent  gas.^^ 

A  further  limiting  factor  in  the  application  of  the  model  presented  in 
Refs.  11  and  12  to  these  experiments  is  that  the  model  does  not  include  the 
mechanisms  vdilch  sustain  the  plasma.  Because  these  mechanisms  are  not 
included,  the  model  cannot  predict  the  voltages  necessaiy  to  initiate  or 
maintain  a  rf  discharge.  These  mechanisms  may  be  important  in 
determining  the  ion  bombardment  energy  distribution.  The  Importance  of 
the  sustaining  mechanisms  is  shown  by  the  changes  in  the  voltage  necessary 
to  maintain  a  rf  discharge  at  different  frequencies. 

Furthermore,  most  of  the  data  presented,  was  collected  at  higher 
frequencies  where  the  assumptions  used  in  the  derivation  of  the  model 
presented  in  Refs.  1 1  and  12  are  no  longer  strictly  valid.  These  assumptions 
also  limit  the  lower  limit  of  the  frequencies  for  which  the  model  is  valid. 
However,  using  the  measurements  that  have  been  performed,  it  is  possible 
to  qualitatively  predict  the  behavior  of  the  ion  bombardment  energy 
distribution  at  pressures  below  100  mTorr  using  the  model.  Furthermore, 
the  measurements  at  higher  pressures  can  be  qualitatively  predicted  by 
modifying  the  ion  bombardment  energy  distributions  predicted  by  the 
model  to  include  the  effect  of  ion-neutral  collisions  in  the  sheath. 


l^see  Ref.  15. 
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The  experiments  that  have  been  performed  have  yielded  information 
which  can  be  used  to  improve  the  model.  These  experiments  have 
confirmed  that  the  model  yields  reasonable  results  for  the  ion  bombardment 
energy  distributions  in  a  rf  neon  discharge.  The  experiments  also  show  areas 
in  which  the  model  needs  to  be  improved. 
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SECTION  V 


A  COLUSIONLESS  PLASBCA  SHEATH  MODEL  FOR  A  PLASBiA  CONTAINING 

NEGATIVE  IONS 
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A  Collisioiiless  Plasma  Sheath  Model  for  a  Plasma  Containing 

Negative  Ions 
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ABSTRACT 

A  sheath  model  Is  derived  for  a  plasma  containing  negative  ions  which 
is  in  contact  with  a  surface.  The  plasma  is  assumed  to  consist  of  a  singly 
charged  positive  ion  species,  a  singly  charged  negative  ion  species, 
electrons,  and  neutrals.  The  model  is  derived  assuming  that  the  entire 
potential  difference  between  the  surface  and  the  bulk  of  the  plasma  occurs 
across  a  sheath  region.  The  thickness  of  the  sheath  region  is  assumed  to  be 
small  with  respect  to  the  mean  free  paths  of  the  particles  Involved,  so  that 
in  the  sheath  region  collisions  and  particle  production  may  be  ignored.  It  is 
shown  that  a  criterion  exists  which  must  be  satisfied  in  order  for  the  model 
to  be  physically  realistic.  This  criterion  is  then  utilized,  under  appropriate 
assumptions  about  the  velocity  distribution  functions  of  positive  ions, 
negative  ions,  and  electrons  entering  the  sheath,  to  derive  expressions  for 
the  conductivity  and  capacitance  of  a  plasma  sheath.  Numerical  results  for 
these  expressions  are  presented  for  cases  which  are  representative  of  a 
typical  plasma  containing  negative  ions.  Implications  of  the  model  for  rf  and 
dc  discharges  are  discussed  including  an  explanation  of  the  confinement  of 
negative  ions  inside  a  plasma  volume. 


I.  INTRODUCTION 

The  study  of  the  properties  of  plasmas  and  plasma  sheaths  has  a  long 
history.  In  order  to  model  the  complete  plasma  phenomena,  the  plasma  has 
historically  been  divided  into  two  regions.  The  first  region  is  the  plasma 
region,  where  the  effect  of  any  space  charge  on  the  electric  potential  may  be 
neglected.  The  second  is  a  sheath  region,  where  the  effect  of  ion  and 
electron  generation  may  be  neglected.  Much  of  the  basis  for  the  modem 
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theory  of  a  pleisma  can  be  found  in  a  paper  by  Tonks  and  Langmuir.^  Their 
paper  was  mainly  concerned  with  the  case  where  the  mean  free  path  of  the 
positive  ions  was  large  with  respect  to  the  dimensions  of  the  container,  so 
that  it  could  be  assumed  that  the  positive  ions  undergo  no  collisions  from 
the  time  they  are  created  until  the  time  they  recombine  at  the  surface.  It 
was  also  assumed  that  the  positive  ions  are  created  with  zero  initial  velocity. 
Using  these  two  assumptions.  Tonks  and  Langmuir  were  able  to  derive  an 
expression  for  the  spatially  dependent  velocity  distribution  and  density  of 
the  positive  ions  as  a  function  of  the  electric  potential  and  the  rate  of 
production  of  the  positive  ions.  By  then  further  assuming  that  the  electrons 
have  a  Maxwellian  velocity  distribution  throughout  the  plasma,  they  derived 
an  equation  (to  subsequently  be  referred  to  as  Tonks  and  Langmuir's  plasma- 
sheath  equation)  for  the  potential  at  any  point  in  the  plasma. 

Tonks  and  Langmuir  were  able  to  approximate  the  solution  to  their 
plasma-sheath  equation  in  the  plasma  region  by  means  of  a  series  expansion 
under  the  assumption  that  the  effect  of  the  space  charge  on  the  electric 
potential  could  be  neglected.  Their  paper  also  included  some  work  for  cases 
in  which  collisions  and  space  charge  effects  could  no  longer  be  neglected. 

The  plasma  sheath  region  was  modelled  by  Bohm.2  He  assumed  that  at 
the  plasma  sheath  edge,  the  velocity  distribution  of  the  positive  ions  is 
monoenergetic  and  that  the  positive  ions  undergo  no  collisions  within  the 
sheath  region.  This  assumption  can  be  valid  in  the  plasma-sheath,  even  for 
plasmas  where  collisions  involving  the  positive  ions  in  the  bulk  of  the  plasma 
cannot  be  Ignored.  This  is  an  idea  which  was  in  part  expressed  earlier  by 
Tonks  and  Langmuir.^  The  assumption  that  the  positive  ions  undergo  no 
collisions  is  valid  if  the  sheath  is  smaller  than  the  mean  free  path  of  the 
positive  ions.  Bohm  also  ignored  any  changes  in  the  potential  within  the 
plasma  region.  Bohm  then  derived  a  criterion  which  the  positive  ions  must 
satisfy  in  order  for  the  model  to  be  physically  realistic.  Bohm's  criterion  was 
that  in  a  plasma  consisting  of  a  singly  charged  positive  ion  species. 


^L.  Tonks  and  I.  Langmuir,  Phys.  Rev.  34.  876  (1929). 

^Davld  Bohm,  In  The  Chamcteristics  of  Electrical  Discharges  in  Magnetic  Fields,  edited  by 
A.  Gutherle,  and  R  K.  Wakerllng  (McGraw-HlU.  New  York.  1949)  chapter  3. 

^sce  Ref.  1.  page  902,  but  there  Is  a  wide  range  of  pressures  for  which  the  Ion  may  drift  in  the 
plasma  yet  faU  freely  through  the  greater  part  of  a  thin  sheath.”  Also,  see  page  901-906  of 
Ref  1  for  a  solution  in  the  sheath  region. 
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electrons,  and  neutrals,  monoenergetic  positive  ions  must  enter  the  sheath 
from  the  plasma  with  a  kinetic  energy  greater  than  or  equal  to  half  the 
electron  thermal  energy. 

Harrison  and  Thompson^  attempted  to  generalize  the  criterion  given 
by  Bohm  to  the  case  of  a  distribution  of  wlocities  for  the  positive  ions. 
However,  as  pointed  out  by  Hall.^  Harrison  and  Thompson's  generalization  is 
not  valid  for  an  arbitrary  distribution  as  was  claimed.  Hall^  also  published 
some  further  important  comments  on  the  misuse  of  the  Bohm  criterion  with 
regard  to  the  concept  of  plasma  stability  and  plasma  oscillations. 

The  concept  of  the  plasma  and  sheath  regions  was  revisited  by  Caruso 
and  Cavallere.'^  Their  contribution  was  to  show  that  the  sheath  region  and 
plasma  region  solutions  are  a  special  form  of  asymptotic  solutions  of  Tonks 
and  Langmuir's  plasma-sheath  equation.  Self^  solved  Tonks  and  Langmuir's 
plasma-sheath  equation  numerically  without  neglecting  the  space  charge 
effects  or  dividing  the  equation  into  sheath  and  plasma  regions.  Selfs 
method  of  solution  eliminates  the  problems  vdiich  arise  when  one  attempts 
to  match  an  asymptotic  plasma  solution  to  an  asymptotic  sheath  solution. 
However,  since  Tonks  and  Langmuir's  plasma  sheath  equation  assumes  that 
the  positive  ions  undergo  no  collisions  after  creation  (i.e.,  a  coUlslonless 
plasma  as  well  as  a  collisionless  sheath).  Selfs  solution  is  limited  in  its 
application  to  actual  dc  discharges. 

All  of  these  previous  studies  Involved  plasma  systems  consisting  of  at 
most  three  particles:  a  single  positive  ion  species,  electrons,  and  neutrals. 
However,  plasmas  containing  negative  ion  species  are  of  current  interest 
due  to  the  use  of  electronegative  gases  in  plasma  processes  used  by  the 
microelectronics  industry.  Negative  ions  are  also  of  current  interest  in  ion 
sources  used  for  negative  ion  beam  production.  In  this  paper,  a  four  particle 
plasma  vdilch  consists  of  a  single  positive  ion  species,  a  sin^e  negative  ion 
species,  electrons,  and  neutrals  will  be  considered.  The  interaction  of  the 
plasma  with  a  surface  which  is  in  contact  with  the  plasma,  leading  to  the 
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formation  of  a  sheath  region,  will  be  modelled.  An  overview  of  the  model  will 
be  presented  In  the  beginning  of  the  next  section.  Subsequently,  the 
assumptions  vdilch  are  made  In  deriving  the  model  and  the  derivation  of  the 
model  will  be  presented.  Finally,  the  applicability  of  the  model's  assumptions 
to  actual  discharges  and  the  implications  of  the  model  for  rf  and  dc 
discharges  will  be  discussed. 

It  will  be  shown  that  even  for  the  situation  where  the  density  of 
negative  ions  in  the  bulk  of  the  plasma  is  much  larger  than  the  density  of 
electrons  in  the  bulk  of  the  plasma,  a  space  charge  sheath  develops  between 
the  plasma  and  a  surface  in  contact  with  the  plasma.  As  a  consequence  of 
this  sheath,  the  surface  still  floats  at  a  significant  negative  potential  with 
respect  to  the  plasma.  Therefore,  the  negative  ions  are  essentially  trapped 
in  the  plasma  by  the  floating  potential  at  the  surface  and,  in  most  instances, 
their  density  in  the  bulk  of  the  plasma  increases  until  volume  losses  equal 
volume  production. 


n.  SHEATH  MODEL  AND  THE  SHEATH  CRITERION 

A  surface  in  contact  with  a  plasma  results  in  a  disturbance  of  the 
plasma.  The  surface  acts  as  a  recombination  center  for  the  positive  ions, 
negative  ions,  and  electrons.  This  results  in  a  flux  of  positive  ions,  negative 
ions,  and  electrons  to  the  surface  and  a  flux  of  neutral  particles  from  the 
surface  back  into  the  plasma.  In  steady  state,  the  loss  of  charged  particles  at 
the  surface  must  be  balanced  by  the  production  of  ions  and  electrons  within 
the  body  of  the  plasma  and  any  net  flux  of  charged  particles  or  electric 
current  collected  at  the  surface  must  be  balanced  by  an  electric  current 
passing  through  an  external  circuit  back  into  the  body  of  the  plasma  via 
another  electrode  surface.  It  is  well  known  that  in  the  case  where  there  are 
no  negative  ions,  if  the  net  current  collected  at  the  surface  Is  zero,  then  the 
surface  will  float  at  a  negative  electric  potential  with  respect  to  the  bulk  of 
the  plasma.  This  potential  is  required  to  reduce  the  value  of  the  electron 
current  collected  at  the  surface  to  the  value  of  the  positive  ion  current 
collected  and  is  due  to  the  differences  in  the  masses  and  average  energies  of 
these  particles.  It  Is  also  expected  that  most  of  the  potential  difference 
between  the  plasma  and  the  surface  will  be  confined  to  a  small  region,  the 
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sheath,  due  to  the  property  that  a  plasma  tends  to  shield  Itself  from  electric 
fields  so  as  to  remain  field  free. 


Plasma  Sheath  Surface 


n+;o.  n€;o,  n.;o  ->x  n+(x),  ne(x).  n-(x)  o 


n+;o.  n€;o,  n.;o  ->x  n+(x),  ne(x).  n-(x)  o 

X  =  0  X  =  Xs 

V(0)  =  0  V(xs)  =  Vs 

E(0)  =  0  Efxs) 

p(0)  *  0  pCxs) 


FIG.  1.  Diagram  of  the  sheath  model  presented  in  section  II. 


The  picture  of  the  plasma  T^ich  will  be  used  in  this  model  is  shown  in 
Fig.  1.  In  this  model,  there  is  an  arbitrary  distinction  between  the  plasma 
and  sheath  regions.  The  plasma  region  supplies  charged  particles  to  balance 
the  charged  particles  which  are  lost  to  the  surface.  The  plasma  region  is 
shielded  from  the  potential  of  the  surface  by  the  sheath  region  £uid  is  a 
relatively  neutral,  field-free,  and  equipotential  region.  However,  in  the 
sheath  region,  the  charge  density,  the  electric  field,  and  electric  potential 
undergo  large  changes.  At  the  surface,  charged  particles  are  neutralized  by 
transferring  their  charge  to  the  surface.  The  surface  will  in  general  have  a 
surface  charge  residing  on  it. 


A.  Assumptions 

In  deriving  this  model,  the  following  assumptions  will  be  made.  These 
assumptions  are  similar  to  the  assumptions  used  by  Bohm^.  with  additional 
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assumptions  regarding  the  negative  ions.  The  applicability  and  validity  of 
these  assumptions  will  be  discussed  in  section  V. 

(1)  The  plasma  consists  of  only  four  types  of  particles:  one  species  of 
singly  charged  positive  ions,  one  species  of  singly  charged  negative  ions, 
electrons,  and  neutrals. 

(2)  The  surface  is  planar  and  both  the  plasma  and  surface  are  infinite 
in  extent,  so  that  the  problem  may  be  treated  as  one-dimensional. 

(3)  All  of  the  voltage  difference  between  the  surface  and  the  plasma  is 
contained  in  the  sheath  region.  The  rest  of  the  plasma  is  considered  to  be 
field  free  and  neutral. 

(4)  The  surface  is  assumed  to  be  at  a  negative  potential  with  respect  to 
the  plasma. 

(5)  The  electrons  are  in  thermodynamic  equilibrium  with  themselves 
and  may  be  characterized  by  an  electron  temperature  Te  which  is  spatially 
Independent.  Also,  it  is  assumed  that  negative  ions  are  in  thermod)mamlc 
equilibrium  with  themselves  and  may  be  characterized  by  a  negative  ion 
temperature  T..  It  is  not  required  that  the  electrons  and  negative  ions  be  in 
equilibrium  with  each  other. 

(6)  The  electron  and  negative  ion  densities  in  the  sheath  are  given  by 
a  Boltzmann  density  distribution. 

(7)  The  positive  ions  undergo  no  collisions  while  traversing  the 
sheath. 

(8)  The  positive  ions  enter  the  sheath  region  with  a  monoenergetic 
drift  velocity.  This  implies  that  at  the  plasma  sheath  edge,  the  ions  have  a 
drift  energy  much  greater  than  their  thermal  energy. 

(9)  Any  charged  particle  which  strikes  the  surface  is  neutralized  at  the 
surface.  No  charged  particles  are  reflected  or  emitted  by  the  surface. 

B.  Derivation  of  the  Sheath  Criterion 

The  electric  field  E(x)  at  any  point  x  in  the  sheath  is  given  by  Guass' 

law 


dE(x)  p(x) 

dx  =  Eo 
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where  p(x)  is  the  charge  density  and  eo  is  the  permittivity  of  free  space.  In 
general,  the  charge  density  will  contain  a  term  for  each  ionized  species.  For 
the  case  being  considered  in  this  paper,  the  charge  density  is  given  by 

p(x)  =  e  [n^(x)  -  ng(x)  -  n.(x)l  .  (2) 


where  e  is  the  magnitude  of  the  charge  carried  by  an  electron  and  n+(x). 
ne(x).  and  n.(x)  are  the  number  densities  of  singly  charged  positive  ions, 
electrons.  £ind  singly  charged  negative  ions  respectively.  The  electric  field  is 
related  to  the  electric  potential  V(x)  by 


E(x)  = 


-  dV(x) 
dx 


(3) 


Multiplying  Eqs  .  (1)  and  (3)  together  yields 
.  dE(x)  -  p(x)  dV(x) 

-  to - — sr- 


(4) 


This  equation  can  then  be  integrated  once  with  respect  to  x  giving 


fE(x)^dx.. 


p(x)  dV(x) 


The  left  side  is  easily  Integrated  to  yield 


2- (E"W  -  £"(0))  =  - 


Assumption  (3)  states  that  the  electric  field  in  the  plasma  is 
insignificant  compared  to  its  value  in  the  sheath,  therefore  we  assume 

E(0)  =  0  .  (7) 

Since  the  electric  potential  due  to  a  distribution  of  free  charges  in  a 
collisionless  region  must  be  a  monotonic  function  of  x.  it  is  po. ;  Ible  to 
express  the  charge  density  in  Eq.  (6)  as  a  function  of  V  rather  than  x 
yielding 


(8) 


V(x) 

e2(x)  =  ^  Jp(V)  dV  . 

The  expression  for  tJie  electric  field  E(x)  in  the  sheath  is  easily  found 
by  noting  that  the  electric  field  s  directed  toward  the  surface  throughout 
the  sheath  region  when  the  surface  is  at  a  negative  electric  potential  with 
respect  to  the  plasma.  Thus, 

E(x)  =  Ve^(x)  !  (9) 

Finally,  since  the  square  of  the  electric  field  must  always  be  positive  as  noted 
by  Bohm2.  the  condition 

V(x) 

—  fp(V)  dV  ^  0  (10) 

Go  0 

must  be  satisfied  for  any  negative  value  of  the  potential  in  the  sheath.  It  is 
this  condition  which  will  be  referred  to  as  the  sheath  criterion. 


C.  Calculation  of  the  Particle  Densities 


Using  assumptions  (4)  and  (5),  the  densities  for  the  negative  particles 
in  the  sheath  region  relative  to  their  values  at  the  sheath  edge  are  given  by 
Boltzmann  density  distributions.  Defining  voltage  equivalents  Ve  and  V.  for 
the  temperatures  Te  and  T.  of  the  electrons  and  negative  ions  respectively 
by 


(11) 


and 


(12) 


where  k  is  Boltzmann's  constant.  The  Boltzmann  density  distributions  may 
be  written  as 
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,  ,  rV(x)>. 

ne(x)  =  n«;o  expf-^J  . 


(13) 


and 


n.(x)  =  n.-o  .  (14) 

Here  ne;o  and  n.;o  are  respectively  the  densities  of  the  electrons  and  negative 
ions  at  the  plasma  sheath  edge,  i.e.,  x  =  0.  By  assumptions  (3),  (5),  and  (6), 
these  are  also  the  values  of  these  densities  throughout  the  plasma  region. 

The  density  of  the  accelerated  particles  (I.e.,  positive  ions)  is 
calculated  by  noting  that  assumptions  (7)  and  (8)  require  the  flux  of  the 
positive  ions  to  be  a  constant  value  throughout  the  sheath.  Therefore, 

n+(x)  v+(x)  =  n+.Q  v+.Q .  (15) 


where  v+(x)  and  v+.q  are  the  velocities  of  the  positive  ions  at  a  point  x  in  the 
sheath  and  at  the  plasma  sheath  edge  (i.e.,  x  =  0)  respectively  and  where  n+.Q 
is  the  positive  ion  density  at  the  sheath  edge. 

Since  the  positive  ions  undergo  no  collisions.  Conservation  of  Energy 
can  be  used  to  relate  the  velocity  of  the  particle  to  the  potential  at  any  given 
point  in  the  sheath,  yielding 


m+  v+2(x) 


2  e 


+  e  V(x)  = 


+  e  V(0) 


(16) 


where  m+  is  the  mass  of  the  positive  ion.  Without  loss  of  generality,  the 
potential  at  the  sheath  edge  can  be  set  equal  to  zero  as  indicated  in  Fig.  1 , 
yielding 


v^(x) 


2  e  V(x) 
m+ 


(17) 


Substituting  this  expression  into  Eq.  (15)  results  in  the  following  expression 
relating  the  positive  ion  density  a^  any  given  point  in  the  sheath  to  the 
electric  potential  at  that  point: 


* 


€ 
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n+(x)  = 


SeVM  • 


(18) 


0^  '^m 


By  defining  the  voltage  equivalent  V^.  of  the  Injection  velocity  as 

2 


"  2  e  ’ 

Ekj.  (18)  may  be  rewritten  as 


(19) 


n+(x)  = 


w 

V. 


(20) 


Using  Eqs.  (2).  (13),  (14)  and  (20),  the  total  net  charge  density  at  any 
location  in  the  sheath  Is  given  by 


p(x) 


/V(x)>, 

(~J 

1  -  n.;o  y  J 

> 

(21) 


D.  Evahiating  the  Sheath  Criterion 


It  is  convenient  to  introduce  the  quantity  Uq  as  the  ratio  of  the 
negative  ion  density  to  the  positive  ion  density  at  the  sheath  edge,  i.e.. 


(22) 


Since  assumption  (3)  requires  charge  neutrality  at  the  sheath  edge. 


^+:0  ■  ■  *^-;0  =  ^  • 


(23) 


Therefore,  the  densities  of  the  negative  particles  at  the  plasma  sheath  edge 
may  be  written  in  terms  of  cxq  and  the  positive  ion  density  at  the  plasma 

sheath  edge, 

-  n+:o  (1  -  Co)  (24) 
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and 


n.;o  =  n+..oao.  (25) 

Consequently.  Eq.  (21)  may  be  re-expressed  in  terms  of  ccq  and  the  positive 
ion  density  at  the  sheath  edge  yielding 

p(x)  =  e  n^..o  — .  -  (1  -  cxo)  exp^^^^^  -  Oq  .  (26) 

Substituting  this  expression  into  Eq.  (8)  and  performing  the  integration 
results  in  the  following  expression  for  the  square  of  the  electric  field  at  any 
point  in  the  sheath; 


In  order  to  satisfy  the  sheath  criterion  of  Eq.  (10),  the  ion  injection 
velocity  is  required  to  be  such  that  the  above  expression  is  always  non- 
negative  for  any  negative  value  of  the  potential.  This  requires  that  the  term 
enclosed  in  braces  in  Eq.  (26)  be  positive  for  any  negative  value  of  the 
potential,  i.e., 

2V,  -^1  ij  +  (1  -  do)  V,[exp^^j-  l] 

+  OoV.^exp^-^^^^J- ij  ^  0  .  (28) 

Solving  inequality  (28)  for  the  square  root  yields 

V^‘^^^2^{(1  -®o)Vc[l  -exp^^^|  +  ctoV.[l  -  exp^^^J 
+  1  .  (29) 
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Since  both  sides  of  E^.  (29)  are  positive  for  all  negative  values  of  the 
potential,  one  can  multiply  the  inequality  by  the  voltage  equivalent  of  the 
injection  velocity  and  then  square  both  sides.  Collecting  the  terms  which 
involve  the  voltage  equivalent  of  the  injection  velocity  results  in  the 
inequality 

{  -  V(x)  -  F[V(x)]}  S  .  (30) 


where 

FIV(x)|  .  (1  -  tto)  V.  [l  Oo V.  [l  -  .  (31) 


The  term  which  multiplies  the  voltage  equivalent  of  the  injection  velocity  on 
the  left  hand  side  of  Eq.  (30)  is  shown  In  Appendix  A  to  be  non-negative  for 
all  negative  values  of  the  potential.  Therefore,  the  voltage  equivalent  of  the 
injection  velocity  for  any  negative  value  of  the  potential  must  satisfy 


F^[V(x)] 

4  {  -  V(x)  -  FlV(x)l]  • 


(32) 


The  minimum  value  of  the  injection  velocity  V+jmm  required  by  Eq.  (31)  is 
equal  to  the  maximum  value  of  the  right  hand  side  of  Eq.  (31)  over  all 
negative  values  for  the  potential. 

Finally,  it  can  be  shown  that  for  the  case  of  no  negative  ion  species 
(l.e.,  (Xo=0),  the  method  used  by  Bohm2  in  deriving  his  classical  sheath 
criterion  utilizing  a  Taylor's  series  expansion  is  equivalent  to  requiring  that 
the  £^.  (32)  hold  at  the  sheath  edge.  In  particular,  since  Eq.  (32)  must  be 
true  for  any  negative  value  of  the  potential.  It  must  be  true  at  the  sheath 
edge  (i.e,,  x  =  0)  where  the  electric  potential  vanishes.  In  this  limit,  as 
shown  in  Appendix  B.  the  inequality  reduces  to 

^ - T’ 

(l-tto)  ^+00^ 


This  Is  a  necessary  condition  for  Ekj.  (10)  to  be  satisfied.  In  addition  for  the 
case  where  ao=0,  as  shown  in  Appendix  G,  Eq.  (33)  is  also  a  sufficient 
condition  for  the  sheath  criterion  of  Ekj.  (10)  to  hold.  Then  the  maximum 
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value  of  the  right  hand  side  of  Eiq.  (32)  does  occur  at  the  plasma  sheath  edge 
and  has  the  value  V+;min=Ve/2  when  ao=0,  which  Is  the  result  obtained  by 
Bohm.  However,  In  the  case  of  negative  Ion  species,  the  mavimiim  value  of 
the  right  hand  side  of  Eq.  (32)  does  not  occur  at  the  plasma  sheath  edge 
where  the  potential  vanishes  and  Eq.  (33).  while  necessary.  Is  not  sufficient 
to  guarantee  that  the  sheath  criterion  given  by  Eq.  (10)  Is  indeed  satisfied. 

m.  SHEATH  CHARACTERISTICS 


A.  Current  Characteristic  and  the  Floating  Potential 

Based  on  the  sheath  model  presented  In  the  previous  sections,  it  is 
possible  to  derive  the  conduction  current  density  through  the  sheath  as  a 
function  of  the  potential  across  the  sheath.  The  conduction  current  density 
due  to  the  positive  Ions  J+  is  independent  of  the  voltage  across  the  sheath, 

and  is  given  by 

fTeTT 

J+*en+.oV+.o*en+:o-y-]5j^.  (34) 


Using  assumptions  (4),  (5),  (6),  and  (9),  the  electron  and  negative  ion 
conduction  currents,  JelV*)  and  J.fVJ  respectively,  are  a  function  of  the 
voltage  across  the  sheath  Vg  and  are  given  by 


Je(Vj  a  -  e  ne;o 


2  7C  m^ 


=  -  e  n+;o  (1 


2  7C  m^ 


(35) 


and 


J-(V8)  a  -  e  n..o 


e^r- 
n  m. 


a  -  e 


f'e  V:-™  fV^\ 


(36) 
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where  mg  and  m.  are  the  electron  and  negative  ion  masses  respectively.  The 
total  conduction  current  density  Jc(Vg)  is  the  sum  of  the  contributions  from 
each  type  of  particle  yielding  the  conduction  current  characteristic 


(37) 


The  floating  potential  Vf  is  then  the  voltage  for  vdilch  the  net 
conduction  current  collected  at  the  surface  is  zero,  l.e.. 


Jc(V^  =  0  . 


(38) 


B.  Capacitance  Chaiacterlstic 

It  is  also  possible  to  calculate  the  capacitance  of  the  sheath.  As  a 
consequence  of  Gauss'  law,  the  surface  charge  density  o  residing  on  the 
surface  bounding  the  plasma  is 

o  =  -  eo  E(Vs)  .  (39) 

If  the  voltage  across  the  sheath  is  time  dependent  then  the  charge  residing 
on  the  surface  Avlll  be  time  dependent,  resulting  in  a  displacement  current 
Jd.  given  by 


_do  dEfVJ  dEfVJdV, 

Jd  =  dt  “  ■  dt  -  ■  ^0  jjVg  dt  * 

From  this  equation,  a  capacitance  per  unit  area  c»  due  to  the  plasma  sheath 
may  be  defined  as 

dEfVJ 

^8  “  ■  6o  (jv  (41) 

8 

The  capacitance  may  be  evaluated  by  noting  that  the  electric  field  at  the 
surface  is  given  by 
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(42) 


Then,  the  resulting  derivative  of  the  electric  field  at  the  surface  with 
respect  to  the  voltage  across  the  sheath  is 


Consequently,  using  Eqs.  (41)  through  (43),  the  capacitance  per  unit  area 
can  be  written  as 


p(Vs) 

=  E(V3)  • 


(44) 


IV.  NUMERICAL  EXAMPLES  OP  SHEATH  CHARACTERISTICS 

In  order  to  perform  numerical  calculations  of  the  sheath 
characteristics,  the  following  values  for  certain  parameters  of  the  plasma 
system  were  assumed: 

m+  =  40  amu, 

m.  =  40  amu, 

n+o  =  0.6  xlO^o  particles/cm3, 

T.  =  500  ‘’K.  and 

Te  =  23,200  “K. 

The  csdculatlons  were  performed  for  various  values  of  oq. 

Figure  2  shows  the  behavior  of  Eq.  (27)  for  E2(V(x))  as  a  ftmctlon  of 
-V(x)  with  (xo*0,3  and  for  several  different  values  of  V+.  In  this  case, 
V+;min=0. 18879  V,  which  is  the  smallest  value  for  which  Eq.  (32)  and  the 
sheath  criterion  of  Eq.  (10)  hold.  As  can  be  seen  from  Fig.  2,  for  values  of 
V+<V+;niin.  there  are  certain  values  of  V(^d  for  vdiich  the  square  of  the 
electric  potential  would  be  negative.  Therefore  choosing  V+<V+;niin  is 
physically  unrealistic.  However,  for  all  values  of  V+>V+;niin  the  square  of  the 


electric  field  is  always  positive  and  hence,  according  to  the  sheath  criterion, 
the  restriction  is  that  must  be  chosen  such  that  V4^V4.;min.  In  the 
following  discussion  of  the  sheath  characteristics,  V^.  will  be  chosen  such 
that  V+=V+;niiu. 


*vw  (V) 

FIG.  2.  The  behavior  of  Eg.  (27)  for  ocQgQ.S  and  for  various  values  of  V+. 


Figure  3  shows  the  dependence  of  V+.min  on  ao  as  determined  by  the 
maximum  value  of  the  right  hand  side  of  Eq.  (32).  For  the  case  of  ao=0.3 
where  V+;inin=0. 18879  V,  this  corresponds  to  an  E2(V)  versus  -VM  curve  (as 
in  Fig.  2)  that  has  the  value  zero  at  -V(x)=0  V,  reaches  a  local  maximum  near 
-V(x)=0.5  V,  decreases  to  zero  near  -V(x)=1.25  V,  and  then  increases  as 
-V(x)  Increases.  Figure  3  also  shows  the  value  of  the  right  hand  side  of 
Eq.  (33)  as  a  function  of  oq.  This  value  requires  that  the  slope  of  the  E2(V) 
versus  -V(x)  curve  at  -V(x)=0  be  nonnegative.  Consider  the  case  for  ao=0.3. 
where  the  value  of  the  right  hand  side  of  Eq.  (33)  equals  0.06840.  For  this 
situation,  the  V+=0.15  curve  in  Fig.  2  is  positive  near  -V(x)=0  since  0.15  is 
greater  than  0.06840.  However,  this  curve  is  negative  for  values  of  -V(x) 
between  approximately  0.75  and  3.5  V.  This  clearly  shows  that  while  the 
criterion  given  in  Eq.  (33)  is  necessary  for  the  solution  to  be  physically 
reallsUc,  It  is  not  always  sufficient.  From  Fig.  3.  one  can  see  that  the 
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criterion  given  in  Eq.  (33)  Is  sufficient  for  cases  in  which  oto  Is  greater  than 
about  0.55  for  the  assumed  discharge  parameters.  Furthermore,  there  are 
values  of  ao  near  zero  for  which  this  criterion  Is  sufficient.  However,  since 
the  criterion  given  in  Eq.  (33)  is  not  sufficient  for  all  ao,  the  criterion  of 
Eq.  (32),  which  is  both  necessary  and  sufficient  for  all  oo.  must  be  used. 


“o 

FIG.  3.  The  dependence  of  V+;min  and  the  right  hand  side  of  Eq.  (33)  on  gp. 


Figures  4a,  4b,  and  4c  show  the  dependence  of  the  electric  potential, 
charge  density,  and  electric  field  respectively  for  ao=0.3  and  V+=V+;niin-  It 
should  be  noted  that  in  order  to  solve  the  equations  for  the  electric  field  and 
electric  potential  numerically,  it  Is  necessary  to  start  at  the  surface  (l.e., 
x=Xs)  where  the  value  of  the  electric  potential  is  Vg  and  the  value  of  the 
electric  field  is  given  by  Eqs.  (27)  and  (9),  and  Integrate  towards  the  plasma. 
This  procedure  Is  required  since,  at  V(x)=0,  all  derivatives  of  V(x)  with 
respect  to  x  vanish.  Thus,  the  electric  potential,  net  charge  density,  and 
electric  field  only  approach  zero  asjonptotically.  This  results  in  the  use  of  x- 
Xg  as  the  spatial  coordinate  in  Fig.  4a  with  Xg  determined  for  Vs=-25  V. 
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FIG.  4a.  The  spatial  distribution  of  the  electric  potential  in  the  sheath  for 
ao=0.3.  V+=V+;niin,  and  V8=-25  V.  In  this  graph,  Xg  is  the  position  of  the 
surface. 


FIG.  4b.  The  dependence  of  the  charge  density  on  the  electric  potential  in 
the  sheath  for  (XpgQ.S  and  V»gV+;min. 
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FIG.  4c.  The  dependence  of  the  electric  field  on  the  electric  potential  In  the 
sheath  for  ao=0.3  and  V^■=V^.;Inln« _ 


Figure  4a  shows  that  most  of  the  potential  In  the  sheath  is  dropped 
across  a  small  region.  For  this  case  with  a  magnitude  of  25  V  across  the 
sheath,  most  of  the  voltage  drop  Is  across  a  region  with  a  dimension  of  about 
1.5  mm.  Figure  4b  shows  the  dependence  of  the  charge  density  on  the 
electric  potential  In  the  sheath.  In  this  example,  there  is  a  region  of  negative 
space  charge  In  the  sheath  near  the  plasma  sheath  boundary  (i.e.,  V(x)=0  V). 
According  to  the  model,  this  may  be  explained  by  noting  that  at  the  plasma 
sheath  boundary,  the  charge  density  is  zero.  After  entering  the  sheath,  the 
negative  ions  are  easily  repelled  by  a  small  negative  electric  potential  and 
this  results  In  the  first  positive  space  charge  region.  As  a  consequence,  the 
positive  ion  density  then  decreases  until  the  positive  Ion  density  becomes 
less  than  the  electron  density  resulting  in  the  negative  space  charge  region. 
Finally,  since  the  electron  density  is  decreasing  exponentially,  the  electron 
density  eventually  becomes  less  than  the  positive  ion  density  resulting  in  a 
positive  space  charge  region  throughout  the  rest  of  the  sheath.  It  should  be 
noted  however,  that  this  detailed  behavior  near  the  plasma  sheath  boundary 
involving  a  negative  space  charge  region  may  not  accurately  reflect  realty 
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because  of  assumption  (3)  used  in  developing  this  model,  as  discussed  in 
Section  V. 

The  current  characteristic  given  by  Eq.  (37)  is  shown  in  Fig.  5a  for 
V+=V+;inin  and  various  values  of  ao-  An  enlargement  of  the  region  near  where 
Jc(Vs)=0  is  shown  in  Fig.  5b.  It  can  be  noted  from  this  figure  that  a  relatively 
large  density  of  negative  ions  with  respect  to  the  density  of  electrons  (i.e.,  ao 
near  1)  is  required  before  there  is  a  significant  effect  on  the  floating 
potential.  Even  for  ao=0.99.  the  floating  potential  is  still  -4.5  V.  However, 
the  introduction  of  even  a  small  density  of  negative  ions  reduces  the  positive 
ion  saturation  current  due  to  the  dependence  of  the  positive  ion  current  on 
V+imin.  which  decreases  rapidly  with  increasing  ao  near  ao=0.  The 
capacitance  characteristic  is  shown  in  Fig.  6.  This  characteristic  is  also 
noticeably  effected  by  changes  in  the  relative  density  of  negative  ions  near 
ao=0.  It  should  also  be  noted  that  the  capacitance  is  a  nonlinear  function  of 
the  voltage  across  the  sheath. 


FIG.  5a.The  conduction  current  density  collected  at  the  surface  as  a  function 
of  the  voltage  between  the  surface  and  the  sheath  for  V+=V+;min  and  for 

different  values  of  oo. 
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FIG.  5b.  Enlargement  of  Fig.  5a.  highlighting  the  region  near  where  Jc(V8)=0. 
The  voltage  V,  for  which  Jc(Va)=0  is  the  value  of  the  floating  potential,  Vf. 
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FIG.  6.  The  capacitance  per  unit  area  of  the  sheath  as  a  function  of  the 
voltage  between  the  surface  and  the  sheath  for  V+=V+;niin  and  for  different 

values  of  oq. 
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V.  APPUCABIUTT  AND  EXTENSIONS  OP  THE  MODEL 


The  assumptions  used  in  the  derivation  of  this  collisionless  plasma 
sheath  model  have  some  important  consequences  with  regard  to  its 
application  and  extensions.  Though  it  has  not  been  explicitly  shown,  it  can 
be  noted  that  it  is  possible  to  allow  any  of  the  ionic  species  to  be  multiply 
charged  by  making  appropriate  changes  in  the  previous  derivation.  In 
addition, while  this  model  can  not  readily  be  extended  to  a  case  involving 
more  than  one  positive  ion  species,  it  can  be  extended  to  a  situation 
involving  more  than  one  negative  ion  species.  Thus,  assumption  (1)  of  the 
model  may  be  relaxed  to  include  any  plasma  consisting  of  neutrals, 
electrons,  single  multiply  charged  positive  ion  species,  and  multiple 
multiply  charged  negative  ion  species. 

Assumption  (2),  which  allows  the  problem  to  be  treated  as  one- 
dimensional,  is  valid  only  if  the  thickness  of  the  sheath  is  much  smaller  than 
the  dimensions  of  the  system.  Assumption  (4)  restricts  the  applicability  of 
the  model  to  cases  where  the  surface  is  at  a  negative  potential  with  respect 
to  the  plasma.  In  addition,  in  order  for  assumptions  (5)  and  (6)  to  be  valid, 
the  magnitude  of  the  electric  potential  at  the  surface  must  be  significantly 
larger  than  the  voltage  equivalent  of  either  the  negative  ion  or  electron 
temperatures.  If  this  condition  is  satisfied,  the  loss  of  electrons  and  negative 
ions  from  the  bulk  of  the  plasma  will  be  si  lall  and  their  velocity  distributions 
should  be  only  slightly  displaced  from  an  equilibrium  situation.  Very  near  the 
surface  however,  there  will  sJll  be  some  departure  from  assumption  (6), 
since  the  Boltzmann  density  distribution  assumes  an  equal  number  of 
particles  moving  both  toward  and  away  from  the  surface,  v/hile,  in  reality, 
near  the  surface  particles  move  predominantly  toward  the  surface.  However, 
since  near  a  surface  at  a  large  negative  potential,  the  electron  and  negative 
ion  densities  are  very  small  with  respect  to  the  positive  ion  density,  any 
deviations  in  their  contributions  to  the  net  charge  density  may  be  neglected. 

Assumption  (7)  is  only  satisfied  if  the  sheath  region  Is  small  with 
respect  to  the  mean  free  path  of  a  positive  ion.  Assumption  (9)  assumes  a 
perfectly  absorbing  and  nonemitting  surface.  It  should  be  noted  that  it  Is 
possible  to  perform  a  derivation  equivalent  to  the  one  presented  here  In 
which  the  surface  Is  assumed  to  have  a  nonzero  reflection  or  emission 
coefficient. 
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Assumptions  (3)  and  (8)  are  the  assumptions  which  most  limit  the 
model.  First,  In  order  for  the  velocity  distribution  of  the  positive  ions  to  be 
considered  monoenergetlc  as  in  assumption  (8),  their  average  drift  energy 
must  be  much  greater  than  their  thermal  energy  In  the  bulk  of  the  plasma. 
This  means  that  V+  must  be  much  greater  than  kT+/e,  where  T+  is  the 
thermal  temperature  of  the  positive  ions.  From  Fig.  3,  it  is  clear  that  when 
Ve  is  much  greater  than  V+,  this  condition  Is  satisfied,  except  as  one 
approaches  the  situation  where  there  are  no  electrons  (i.e.  ao=l)  and  T.  is 
less  than  or  equal  to  T+.  However  in  reality,  even  if  the  positive  ions  enter 
the  sheath  with  an  average  drift  energy  which  is  much  larger  than  the 
thermal  energy  of  these  ions  in  the  bulk  of  the  plasma,  the  width  of  the 
velocity  distribution  function  of  these  ions  may  be  much  broader  than  is 
expected  from  the  thermal  energy  of  the  ions  in  the  bulk  of  the  plasma.  This 
is  result  of  the  fact  that  the  ions  are  flowing  from  the  point  at  which  they  are 
created  to  the  plasma  sheath  boundary  and  any  collisions  that  these  ions 
undergo  before  they  reach  the  plasma  sheath  boundary  will  broaden  their 
energy  distribution.  However,  this  assumption  will  still  be  valid  throughout 
most  of  the  sheath  region,  except  near  the  plasma  sheath  boundary,  since  in 
the  sheath  the  ions  gain  large  amounts  of  energy  relative  to  the  kinetic 
energy  with  which  the  ions  are  injected  into  the  sheath. 

The  above  discussion  points  out  the  fact  that  this  model  does  not 
address  the  origin  of  the  drift  velocity  other  than  noting  that  it  is  due  to  the 
fact  th^t  positive  ions  are  created  in  the  plasma  and  are  lost  at  the  surface. 
Other  authors  have  attempted  to  model  the  origin  of  the  drift  velocity 
through  the  use  of  a  presheath  concept.  However,  as  implied  by 
assumption  (3),  the  use  of  a  presheath  has  been  avoided  in  this  discussion.  It 
should  be  noted  that  by  using  the  equation  for  the  electric  potential  in  the 
plasma  sheath  resulting  from  this  model  under  the  boundary  conditions 
Indicated  in  Fig.  1,  one  finds  that  the  solution  for  the  el''ctric  potential  only 
asymptotically  approaches  zero  as  one  proceeds  towards  the  plasma  from 
the  surface.  Hence,  the  sheath  would  be  infinite  in  extent  if  it  is  required 
that  the  electric  field  and  the  net  charge  density  be  zero  at  the  plasma 
sheath  boundary,  as  was  assumed  in  this  derivation.  However,  in  reality,  the 
electric  field  and  the  net  charge  density  have  a  small  but  nonzero  value  at 
the  plasma  sheath  edge  and  the  resulting  solution  of  the  sheath  then  has  a 


82 


finite  size.  In  either  case,  however,  most  of  the  potential  in  the  sheath  is 
dropped  over  a  finite  distance. 

VI.  MODEL  IMPUCATIONS  FOR  DC  AND  RF  DISCHARGES 

This  model  has  several  Important  implications  for  dc  and  rf  discharges 
in  electronegative  gases.  For  dc  discharges,  the  model  predicts  that  in  the 
positive  column  it  is  possible  to  achieve  a  very  large  density  of  negative  ions. 
As  previously  shown,  the  floating  potential  decreases  slowly  as  the  ratio  of 
the  negative  ion  density  to  the  positive  ion  density  increases.  Hence,  the 
voltage  equivalent  of  the  negative  ion  temperature  is  much  smaller  in 
magnitude  than  the  floating  potential,  even  when  the  negative  ion  density  is 
greater  than  99%  of  the  total  density  of  negatively  charged  particles. 
Therefore,  practically  no  negative  ions  are  able  to  escape  to  the  walls  of  the 
dc  discharge.  The  negative  ions  are  essentially  trapped  in  the  plasma  by  the 
floating  potential  at  the  surface.  Hence,  in  most  cases  involving  negative 
ions,  their  density  in  the  bulk  of  the  plasma  increases  until  volume  losses 
equal  volume  production. 

For  rf  discharges,  it  can  also  be  shown  that  the  negative  ions  are 
trapped  within  the  plasma  in  a  similar  manner.  By  the  same  argument  as 
given  above,  practically  no  negative  ions  may  escape  to  any  side  walls  of  a  rf 
reactor  since  these  walls  are  at  the  floating  potential.  In  addition,  using  a 
model  of  a  rf  plasma  that  has  been  developed  by  Metze,  Ernie,  and  Oskam^, 
the  loss  of  negative  ions  to  the  electrodes  may  be  predicted.  This  model 
shows  that  the  magnitude  of  the  electric  potential  across  the  plasma  sheath 
at  an  electrode  is  never  significantly  less  than  magnitude  of  the  floating 
potential.  This  result  is  a  consequence  of  che  current  characteristic  and 
capacitance  characteristic  of  the  plasma  sheath.  When  the  magnitude  of  the 
sheath  voltage  is  only  slightly  less  than  magnitude  of  the  floating  potential, 
very  large  electron  currents  can  be  collected.  However,  when  the  magnitude 
of  the  sheath  potential  is  greater  than  the  magnitude  of  the  floating 
potential,  the  maximum  current  that  can  be  collected  is  the  much  smaller 
positive  ion  saturation  current.  Since  the  net  current  through  the  sheath 
during  an  rf  cycle  is  zero  if  a  blocking  capacitor  is  In  series  with  the 

9a.  Metze,  D.  W.  Ernie,  and  H.  J.  Oskam.  J.  y^pl.  Phys.  60.  3081  (1986). 


electrode,  the  magnitude  of  the  voltage  across  the  sheath  during  that 
portion  of  the  cycle  during  which  a  large  electron  current  can  flow  through 
the  sheath  must  be  only  slightly  smaller  than  the  magnitude  of  the  floating 
potential.  Again,  for  this  argument,  since  the  magnitude  of  the  electrode 
potential  Is  always  near  or  greater  than  the  magnitude  of  the  floating 
potential,  practically  no  negative  Ions  ever  escape  to  the  electrodes.  Hence, 
in  most  rf  plasmas,  the  negative  Ion  density  v;dll  Increase  imtll  volume  losses 
equal  volume  production.  Therefore,  a  very  large  density  of  negative  ions 
may  be  expected  for  plasmas  In  electronegative  gases. 
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APPENDIX  A 

In  this  Appendix,  the  coefficient  of  V+  on  the  left  hand  side  of  Eq.  (30) 
is  shown  to  be  greater  than  or  equal  to  zero,  l.e., 

-  V(x)  -  Frv(x)]  ^  0  .  (Al) 

for  V(x)  ^  0. 

It  will  first  be  useful  to  prove  that  the  function  g(z)  defined  by 
g(z)  =  z  +  exp(  -  z)  -  1  (A2) 

is  non-negative  for  z  non-negative.  To  prove  this  statement,  first  note  that 
g(0)  =  0  .  (A3) 

The  derivative  of  the  function  g(z)  with  respect  to  z  is 

=  1  -  exp(  -  z)  .  (A4) 

Since  z  is  non-negative,  the  derivative  Is  also  non-negative.  Therefore, 
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g(z)  ^  0 


(A5) 


if  z  is  non-negative,  as  desired. 

The  substitution  of  -V(x)  A^e  for  2  hi  Eqs.  (A2)  and  (A5)  is  valid  since,  if 
V(x)  ^  0,  then  -V(xJ/Ve  is  non-negative.  Multiplying  both  sides  of  the 
resulting  inequality  by  Vg  yields 

-  V(x)  -  Ve[l  -  exp^^)]  ^  0  •  (A6) 

Similarly,  using  V.  for  the  negative  ions, 

-  V(x)  -  V.  [  1  -  exp^^^jj  ^  0  .  (A7) 

Multiplying  Eq.  (A6)  by  (1  -  ao).  and  Eq.  (A7)  by  Oo.  and  then  adding,  yields 

-  V(x)  -  (1  -  tto)  Ve  [^1  -  exp(^^^^  -  aoV.|^l  -  exp^'^^j]  >  0  .  (A8) 

But  from  the  definition  of  F(x)  in  Eq.  (31).  Eq.  (A8)  is  equivalent  to  Eq.  (Al). 
Thus,  Eq.  (Al)  is  shown  to  hold  for  all  V(x)  ^  0. 

APPENDIX  B 

In  this  Appendix,  it  will  be  shown  that  Eq.  (33)  is  the  limit  of  Eq.  (32) 
as  the  potential  V(x)  approaches  zero. 

When  V(x)  is  equal  to  zero.  Eq.  (32)  Is  an  indeterminate  form  which 
may  be  evaluated  by  expan'  .ing  the  numerator  and  denominator  in  a  Taylor's 
series  about  the  point  V(x)  =  0.  Expanding  F[V(x)l  yields 

F[V(x))  =  (1  -  tto)  Ve[l-(^1  + 

+  aoV.[l-(^l  +  (Bl) 

where  the  neglected  terms  are  of  a  higher  order  in  V(x).  This  expression 
simplifies  to 
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F[V(x)J  =  -V(x)  -  (1  -  ao) 


V^(x)  ^  V^(x)  , 
2  Vg  "  2  V. 


(B2) 


Substituting  this  expression  Into  Ekj.  (32)  yields 

.  V^(x)  V^(x) 

j^“V(x)  -  (1  ■  oto)  2  V  "  ®o  2  V 

V+  ^  r  v^fx)  V^bd  tT 

4  |-V(x)  -  [-V(x)  -  (1  -  ao)-2^-  Oo  2^+..-J 


(B3) 


Simplifying  this  expression  yields 
_ ^(x)  +  ... 


(1  -  ao)^+ao—  V2(x)  + 


(B4) 


In  the  limit  as  V(x)  approaches  zero,  the  terms  of  higher  order  in  V(x) 
in  the  numerator  and  denominator  may  be  neglected.  Thus,  cancelling  the 
V2(x)  term  from  both  the  numerator  and  denominator  yields  the  result  given 
in  Eq.  (33). 


APPENDIX  C 

In  this  Appendix,  it  will  be  shown  that  for  the  case  of  a  only  one  singly 
charged  positive  ion  species  and  electrons,  Eq.  (33)  Is  a  sufficient  condition 
for  Eq.  (10)  to  be  satlsfled  for  V(x)^0. 

It  will  be  useful  to  first  prove  that  the  function  g(z)  defined  by 

g(z)  =  2  (z  +  exp(  -  z)  -  1]  -  (1  -  exp(-  z)]2  ,  (Cl) 

Is  non-negative  for  z  non-negative.  To  prove  this  statement,  first  note  that 

g(0)  =  0  .  (C2) 

The  derivative  of  the  function  g(z)  with  respect  z  is 

=2(1-  exp(  -  z)J  -  2  (1  -  exp(  -  z)J  exp(  -  z)  ,  (C3) 

2ind  It  may  be  rewritten  as 
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(C4) 


=  2  [1  -  exp{  -  z)]2  . 


This  derivative  Is  clearly  non-negative.  Therefore. 


g(z)  ^  0 


(C5) 


If  z  Is  non-negative,  as  desired. 

Using  Eq.  (Cl),  Eq.  (C5)  may  be  rewritten  as 

2  [z  +  exp(  -  z)  -  1)  ^  [1  -  exp(  -  z)]2  .  (C6) 


for  z  non-negative.  The  left  hand  side  of  Eq.  (C6)  Is  positive  and  the 
Inequality  may  be  rewritten  as 


1  ^  [1  -  exp(  -  z)l2 

^  ^  2  [z  +  exp(  -  z)  -  1]  • 


(C7) 


where  z  Is  non-negative.  Caution  Is  needed  In  performing  this  division  since 
the  numerator  and  denominator  both  go  to  zero  as  z  goes  to  zero.  However, 
the  division  Is  valid  In  the  limit  as  z  goes  to  zero  and  the  value  of  the  right 
hand  side  at  z  =  0  will  be  defined  to  be  this  limit. 

The  substitution  of  -V(x)/Ve  for  z  In  Eq.  (C7)  is  valid  since  If  V(x)  ^  0 
then  -V(x)/Ve  Is  non-negative.  Multiply  both  sides  of  the  resulting  Inequality 
by  Ve/2  yields 


1  -  expj 

rv(x)>, 

r 

4- 

[-  V(x) 

i-V, 

1 

L 

-  expj^ 

V(x)>.-Jl 

V.  Jjl 

(C8) 


The  right  hand  side  Is  clearly  the  same  as  the  right  hand  side  of  Eq.  (32)  for 
the  case  vdien  there  are  no  negative  Ions,  l.e.,  Uq  =  0.  This  proves  that 

Eq.  (33)  Is  a  sufficient  condition  for  the  sheath  criterion  given  by  Eq.  (10)  to 
be  satisfied  for  the  case  when  there  is  only  one  singly  charged  positive  Ion 
species  and  electrons  (l.e.,  (Xo=0). 
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SECTION  VI 


THE  RADIAL  DISTRIBUTIONS  OP  CHARGED  PARTICLE  DENSITIES  AND 
ELECTRIC  FIELD  STRENGTH  IN  THE  POSITIVE  COLUMN 
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Department  of  Electrical  Engineering 
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Minneapolis,  MN  55455 


ABSTRACT 

A  model  is  presented  for  the  radial  distribution  of  ion  and  electron 
densities  and  electric  field  strength  in  the  positive  column  of  a  dc  discharge 
for  a  plasma  consisting  of  a  singly  charged  positive  ion  species,  electrons, 
and  neutrals.  The  set  of  equations  involved  consists  of  the  particle  and 
momentum  conservation  equations  for  the  ions  and  electrons  and  Poisson's 
equation.  Utilizing  this  single  set  of  equations  and  appropriate  assumptions, 
this  model  has  been  solved,  through  suitable  numerical  techniques,  for 
various  tube  radii,  R,  and  gas  pressures,  p©.  These  calculations  show  the 
development  of  both  the  ambipolar  electric  field  in  the  "bulk"  of  the  positive 
column  and  the  sheath  field  "near"  the  discharge  wall.  The  results  also 
demonstrate  the  existence  of  a  nonzero  difference  between  the  ion  and 
electron  densities  at  the  discharge  axis,  with  an  Increase  in  this  difference 
for  decreasing  Rpo. 


I.  INTRODUCTION 

Several  authors  have  reported  theories  concerning  the  radial  distributions  of 
charged  particle  densities  and  the  electric  field  strength  in  the  cylindrical 
positive  column  of  a  direct  current  discharge.  These  theories  were  all 
confined  to  plasmas  containing  a  singly  charged  positive  ion  species  and 
electrons,  with  the  authors  making  various  assumptions  in  order  to  arrive  at 
on  analytic  solution  or,  more  recently,  a  numerical  analysis.  Schottkyi 
developed  an  ambipolar  diffusion  theory  in  1924  in  which  he  assumed 


Present  address:  Honeywell  Systems  and  Research  Center,  3660  Technology  Drive, 
Minneapolis,  MN  55418. 

Iw.  Schottky,  Z.  Phys.  28,  625  (1924). 
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quasineutrallty  of  the  plasma  and  neglected  the  effect  of  charged  particle 
Inertia  on  the  radial  structure  of  the  positive  column.  Subsequent  authors 
have  also  assumed  quasineutrallty  of  the  plasma,  but  included  the  charged 
particle  inertia  term  in  the  conservation  of  momentum  equatlons.2.3  Several 
other  authors  included  Poisson’s  equation  for  calculating  the  radial  electric 
field  strength  inside  the  plasma.^*®  However,  in  order  to  arrive  at  a  tractable 
solution,  they  either  omitted  the  charged  particle  inertia  term^-e  or  the 
effect  of  ion  density  diffusion  on  the  ion  current  density  towards  the 
discharge  walls'^-®.  All  the  authors  have  assumed  a  Maxwellian  distribution  of 
the  random  energy  of  the  ions  and  electrons  with  constant  ion  and  electron 
temperatures. 

This  paper  deals  with  calculations  of  the  radial  charged  particle  density 
distributions,  the  radial  charged  particle  current  densities,  and  the  radial 
electric  field  strength  in  the  cylindrical  positive  column  of  a  dc  discharge. 
The  set  of  equations  involved  consists  of  the  particle  conservation  equations, 
the  momentum  conservation  equations  Including  the  diffusion  and  inertia 
terms,  and  Poisson’s  equation.  The  common  assumption  that  the  energy 
distribution  of  the  random  energies  of  the  electrons  and  ions  is  Maxwellian 
with  temperatures  independent  of  radial  position  will  be  made.  In  the  next 
section  the  the  assumptions  used  in  deriving  the  model  and  the  resulting 
set  of  equations  will  be  given,  along  with  the  method  of  solution  used  in  the 
numerical  calculations.  The  results  obtained  for  realistic  discharge 
parameters  £ind  for  various  values  of  the  discharge  tube  radius,  R,  and  gas 
pressure,  po.  will  be  presented  and  discussed  in  Sec.  III.  The  validity  of  the 
assumptions  used  will  be  compared  with  the  calculated  results  in  order  to 
determine  the  self-consistency  of  the  theory  presented. 


2s.  A.  Self  and  H.  N.  Ewald.  Phys.  Fluids  9.  2486  (1966). 

3h.  W’.  Frtedman.  Phys.  Fluids  10,  2053  (1967). 

P.  Allis  and  D.  J.  Rose.  Phys.  Rev.  93.  84  (1954). 

M.  Cohen.  Phys.  nulds  6,  1492  (1963). 

®I.  M.  Cohen  and  M.  D.  Kruskal.  Phys.  Hulds  8.  920  (1965). 
■^J.  R  Forrest  and  R  N.  Franklin.  J.  Phys.  D  1.  1357  (1968). 
®J.  H.  Ingold,  Phys.  Rulds  15.  75  (1972). 
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n.  THEORETICAL  MODEL 


In  the  first  part  of  this  section,  the  assumptions  used  in  deriving  the 
model  will  be  discussed.  The  resulting  equations  used  for  the  calculations  of 
the  radial  distributions  of  the  charged  particle  densities  and  the  electric 
field  strength  will  then  be  presented.  Finally,  the  botmdaiy  conditions  for 
the  model  and  the  method  of  solution  of  the  set  of  equations  will  be  given  in 
the  last  two  parts  of  this  section. 


A.  Assumptions 

The  model  and  calculations  presented  in  this  paper  refer  to  a  positive 
column  of  a  dc  gaseous  discharge  under  cylindrical  S)rmmetry.  The  plasma 
is  assumed  to  consist  of  a  singly  charged  positive  ion  species,  electrons,  and 
neutrals.  The  following  assumptions  will  be  made  for  the  present  model: 

(1)  The  energy  distribution  of  the  random  energy  of  the  ions  and 
electrons  is  Maxwellian  with  temperatures  Ti  and  Te,  respectively.  This  is 
an  assumption  made  by  previous  authors  and  is  valid,  except  at  locations 
close  to  the  discharge  tube  'vall.  The  latter  point  will  be  discussed  in  more 
detail  when  presenting  the  boundary  conditions  of  the  model. 

(2)  All  quantities  are  independent  of  the  axial  coordinate  z  and  the 
azimuthal  angle  <|>. 

(3)  The  ion  and  electron  temperatures  Ti  and  Te  are  independent  of  the 
radial  coordinate  r.  This  sissumptlon  was  also  made  by  all  previous  authors. 

(4)  The  ionization  frequency  vi  of  the  electrons  is  very  small  with 
respect  to  the  collision  frequencies  for  momentum  transfer  Vmi  and  Vme  of 
the  ions  and  electrons,  respectively.  Therefore,  the  effect  of  ionizing 
collisions  on  the  momentum  loss  (gain)  of  the  charged  particles  can  be 
neglected. 

(5)  The  ion  and  electron  production  is  a  linear  function  of  the  electron 
density,  while  the  only  charged  particle  loss  from  the  plasma  is  by  diffusion 
towards  the  discharge  tube  wall.  The  calculations  can  be  extended  to 
nonlinear  production  (loss)  processes.  However,  the  main  features  of  the 
results  presented  would  not  change  significantly. 
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(6)  The  friction  force  between  the  Ions  and  electrons  can  be  neglected 
with  respect  to  the  friction  force  between  the  charged  particles  and  the 
neutral  particles.  This  is  a  valid  assumption  for  low  degrees  of  ionization. 

(7)  The  positive  column  is  in  steady  state  with  no  net  current  flowing  to 
the  discharge  tube  walls,  so  that  the  radial  ion  and  electron  current 
densities  are  equal. 


B.  Equations 

The  assumptions  given  above  result  In  the  following  set  of  equations: 
1  d 


7^[r  Trlr)]  =  Vi  ne(r) 


(1) 


mi  d 


r  dr 


me  d 


rr^(r) 


“nTfrTJ  ®  ^r(r)  Vmi 


(2) 


r  dr 


r  r^(r) 


L  ne(r)  J 


+  kTe^ne(r)  +  e  ne(r)  Er{r)  =  -  me  Trlr)  Vme  (3) 


f  ^  ^  [ni(r)  -  ne(r)] .  (4) 

Here,  Tflr)  Is  the  radial  particle  flow  density  of  the  Ions  and  electrons:  ni(r) 
and  ne(r)  are  the  Ion  and  electron  particle  densities:  mi  and  me  are  the  Ion 
and  electron  masses;  Eq  Is  the  permittivity  of  free  space:  and  E^(r)  Is  the 
radial  component  of  the  electric  field  strength. 

Equation  (1)  Is  the  continuity  equation  of  the  charged  particles,  while 
Eqs.  (2)  and  (3)  are  the  conservation  of  momentum  equations  (volume  force 
equations)  for  the  ions  and  electrons,  respectively.  Equation  (4)  is  Poisson's 
equation  and  Is  needed  for  a  self-consistent  calculation  of  the  radial  electric 
field  strength.  The  assumption  that  the  ion  and  electron  temperatures  are 
Independent  of  coordinates  eliminates  the  need  to  use  the  conservation  of 
energy  equations  for  the  Ions  and  the  electrons. 
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The  term  on  the  right  hand  side  of  Eqs.  (2)  and  (3)  is  the  volume  friction 
force.  Since  the  radial  drift  velocity  of  the  neutral  particles  is  very  small 
with  respect  o  the  radial  drift  velocity  of  the  ions  and  electrons,  it  is 
neglected  in  this  term.  The  terms  on  the  left  hand  side  of  Kqs.  (2)  and  (3) 
are  the  volume  inertia  force,  the  volume  diffusion  force,  and  the  volume 
electric  field  force,  respectively.  Previous  theories  of  the  radial  properties 
of  the  positive  column  which  incorporated  Poisson’s  equation  have  used 
various  approximations  concerning  the  effect  of  the  different  volume  force 
terms  on  the  radial  particle  density  distributions,  etc.  Several  authors 
neglected  one  or  both  volume  inertia  force  terms, ‘**7  while  one  author 
neglected  the  volume  diffusion  force  for  the  ions®.  The  motivation  for  these 
approximations  result  from  the  fact  that  when  all  three  volume  force  terms 
related  to  the  ions  are  included  in  the  theoretical  model,  an  Instability 
occurs  in  the  numerical  calculations  (for  certain  numerical  techniques) 
when  the  radial  ion  drift  energy  is  equal  to  one  half  the  kinetic  ion  pressure, 
i.e.,  the  radial  ion  drift  velocity  v^i  is  equal  to  the  Isothermal  ion  sound 
speed  Vri=(kTi/mi)i/2,  This  instability  problem  in  the  numerical  techniques 
has  been  discussed  in  the  literature.®*®  The  method  described  below  for 
solving  the  set  of  Eqs.  (l)-(4)  circumvents  this  difficulty. 


C.  Boundary  Conditions 

The  four  functions  to  be  determined  from  Eqs.  (l)-(4)  are  ni(r),  ne(r), 
rr(r),  and  Er(r).  Appropriate  boundary  conditions  are  thus  required  for 
these  functions.  Because  of  symmetry,  both  Frfr)  and  Er{r)  must  be  zero  at 
the  axis  of  the  positive  column,  i.e..  rr(r=0)=0  and  Er(r=0)=0.  The  two 
other  required  boundary  conditions  are  the  values  of  ni(r=0)  and  ne(r=0). 
For  the  model  presented  here,  these  values  will  be  different,  i.e.,  ni(r)?tne(r). 
since  It  will  not  be  assumed  that  the  discharge  plasma  Is  quasineutral.  It 
should  be  noted  that  the  values  of  ni(r=0)  and  ne(r=0)  then  determine  the 
spatial  derivative  of  Er(r)  at  r=0.  Moreover,  the  cylindrical  symmetry  of  the 
positive  column  requires  that  the  radial  derivatives  of  ni(r)  and  ne(r)  equal 
zero  at  r=0. 


®H.  W.  FYledman  and  E.  Levi.  Phys.  Fluids  10.  1499  (19671. 
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The  boundary  condition  for  the  electron  density  at  the  wall  of  the 
discharge  tube  has  been  discussed  by  several  authors.  A  s-ommaiy  of  this 
discussion  is  given  in  Ingold.®  If  it  is  assumed  that  the  velocity  distribution 
of  the  electrons  at  the  wall,  which  is  assumed  to  be  a  perfect  sink  for  the 
electrons,  is  Maxwellian  in  the  forward  direction  and  zero  in  the  backward 
direction,  then  the  following  boundary  condition  must  hold  at  the  wall: 


ne(R)  =  rr(R) 

me  ) 


(5) 


This  condition  is  also  the  boundary  condition  resulting  from  collisionless 
space  charge  theory  and  is  equal  or  very  close  to  the  boundary  condition  of 
other  theories. 


D.  Method  of  Solution 

The  method  of  solution  used  in  this  paper  consists  of  expanding  the 
functions  ni(r).  ne(r),  TfCr),  and  Er(r)  in  power  series  of  the  radial 
coordinate  r,  around  rsO.  The  calculation  of  the  coefficients  of  the  four 
power  series  is  performed  in  the  following  steps: 

(1)  The  boundary  conditions  rr(r=0)=0  and  E^(r=0)=0  as  well  as  the 
symmetry  conditions  (i.e.,  the  radial  derivatives  of  the  charged  particle 
densities  equal  zero  at  r=0)  are  used.  The  values  of  ni(r=0),  ne(r=0),  and  the 
constants  vi,  Vmi,  Vme.  Ti,  and  Te  are  chosen  and  substituted  into  Eqs.  (Il¬ 
ls).  These  values  Avlll  depend  on  the  parameters  of  the  experimental 
situation  being  modeled,  such  as  the  electric  current  density  inside  the 
positive  column,  the  type  of  gas,  the  gas  pressure,  the  tube  radius,  etc,  and 
can  be  appropriately  chosen.  For  example,  the  collision  frequencies  and 
charged  particle  temperatures  may  be  available  from  experimental  data 
presented  in  the  literature  or  can  be  theoretically  calculated. 

(2)  The  coefficients  of  the  power  series  are  determined,  either  by  a 
standard  computer  technique  or,  if  possible,  by  deriving  a  recurrence 
relation  between  the  coefficients.  During  the  calculations,  non-physical 
solutions  for  the  unknown  functions  may  occur.  As  a  result  one  of  the 
parameters,  vi,  is  changed  until  a  physical  solution  having  sufficient  accuracy 
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is  obtained.  The  accuracy  can  be  determined  by  extending  the  power  series 
to  higher  powers  of  r  and  comparing  the  results. 

(3)  The  radius  of  the  discharge  tube  walls  is  determined  from  Eq.  (5).  If 
the  calculated  tube  radius  is  different  from  the  desired  tube  radius,  the 
calculations  are  repeated  using  a  different  value  for  the  net  charged  particle 
density  at  the  axis  An(r=0)=ni(r=0)-ne(r=0).  The  calculations  are  repeated 
until  the  proper  consistency  Is  reached  between  the  tube  radius  for  the 
experimental  conditions  being  modeled  and  the  theoretical  results. 

It  should  be  noted  that  the  set  of  Eqs.  (l)-(4)  can  easily  be  transformed 
into  a  set  of  dimensionless  equations,  with  the  same  method  of  solution 
being  used  to  obtain  the  dependence  of  positive  column  properties  on 
discharge  parameters.  Most  previous  authors  have  used  the  dimensionless 
approach  when  developing  theories  concerning  the  positive  column. 
However,  since  the  present  paper  attempts  to  model  the  properties  of 
positive  column  discharges  for  which  experimental  information  is  available, 
the  dimensionless  approach  has  been  avoided. 


m.  NUBflSRICAL  CALCULATIONS  AND  DISCUSSION 

Numerical  calculations  were  performed  for  the  model  presented  in  Sec. 
II,  utilizing  the  method  of  solution  outlined.  These  calculations  were  for  a 
discharge  in  helium  with  the  following  assumed  parameters: 
mi  =  4  amu  (for  He+), 

Vmi/Po  =  3.0  X  107  s*i  (for  He+), 

Vme/Po  =  2.3  X  109  S'l. 

Ti  =  300  OR, 

ni(r=0)  =  1.0  X  1016  m-l, 

where  the  values  for  Vmi  and  Vme  are  from  Ref.  10.  Results  are  presented 
below  for  different  values  of  R  and  p©.  For  each  selection  of  R  and  po,  the 
appropriate  value  for  T®  was  chosen  from  the  experimental  work  of  Leiby 
and  Rogersii  and  the  values  of  vi  and  ne(r=0)  were  adjusted,  as  indicated  in 


l^J.  H.  Ingold  and  H.  J.  Oskam,  Phys.  Fluids  27. 214  (1984). 

1  ic.  C.  Leil^.  Jr.  and  C.  W.  Rogers,  presciiled  at  the  20th  Gaseous  Electronics  Conference,  Palo 
Alto.  Caltfomla.1967  (unpublished). 
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Sec.  II.  to  give  a  physical  solution  consistent  with  the  chosen  values  of  R  and 
Po- 

Figures  1-3  show  the  results  of  calculations  for  a  discharge  tube  radius 
R=1.0  cm  and  a  gas  pressure  po=5  Torr.  leading  to  the  product  Rpo=5.0  cm- 
Torr.  From  the  results  of  Lelby  and  Rogers,  these  conditions  imply  an 
electron  temperature  Te=26,000  oR.  The  required  net  charged  particle 
density  on  the  axis  is  then  An(r=0)=!6.8  X  10^2  ni-3,  with  an  ionization 
collision  frequency  vi=1.95  X  10^  s-i.  Figure  1  gives  the  plots  of  the  Ion 
density  ni(r),  the  electron  density  ne(r).  and  the  net  charged  particle 
density  An(r).  The  curves  for  ni(r)  and  nc(r)  are  close  to  the  shape  of  the 
zero  order  Bessel  function  Jo(2.4r/R)  resulting  from  the  ambipolar  diffusion 
theory  for  the  positive  column.  ^2  However,  there  is  a  deviation  from  this 
shape,  especially  for  ni(r),  near  the  tube  wall.  This  deviation  Is  Indicative  of 
the  development  of  a  thin  sheath  region  near  the  wall,  which  Is  not  modeled 
by  the  ambipolar  theory.  The  development  of  this  sheath  region  is  also 
Indicated  in  the  curve  of  An(r),  vdilch  Increases  rapidly  near  the  wall.  From 
these  results  it  can  be  seen  that,  while  this  model  predicts  the  development 
of  a  sheath  region,  there  is  no  clear  distinction  between  the  sheath  region 
and  the  bulk  plasma,  as  can  be  expected.  Therefore,  the  assignment  of 
different  regions  of  the  discharge  to  these  two  categories  is  somewhat 
arbitrary,  as  discussed  previously  by  Self  and  Ewald.2.  Finally,  the  instability 
point  discussed  in  Sec.  II  at  which  the  the  radial  ion  drift  velocity  vdi  Is 
equal  to  the  isothermal  Ion  sound  speed  Vji  Is  also  shown  in  Fig.  1 .  This 
point  is  very  near  the  wall  under  these  discharge  conditions. 

The  radial  electric  field  Er(r)  and  the  resulting  electric  potential  V(r)  are 
shown  in  Fig.  2.  For  comparison,  the  ambipolar  radial  electric  field  Ea(r) 
calculated  from  the  ambipolar  diffusion  theory  is  also  plotted  in  this  figure. 
The  value  of  E^(r)  is  relatively  small  near  the  axis,  but  increases  to  larger 
values  near  the  wall.  This  is  again  indicative  of  the  development  of  the 
sheath  near  the  wall.  A  comparison  between  E^fr)  and  Egfr),  indicates  that 
the  two  values  are  equal  near  the  axis,  but  deviate  from  each  other  near  the 
wall.  This  deviation  Is  due  to  the  fact  that  the  ambipolar  theory  does  not 
model  the  sheath  region,  with  Ea(r)  ->+oo  for  r  R  .  The  curve  of  V(r) 
also  shows  the  smooth  transition  from  the  bulk  plasma  region  near  the  axis 


12h.  J.  Oskam,  Philips  Res.  Rep.  13. 335  (1958). 
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Fig.  2.  Numerical  c2ilculations  of  the  radial  electric  field  strength  Er(r)  and 
the  resulting  electric  potential  V(r)  for  Rpo*5.0  cm-Torr.  The  ambipolar 
radial  electric  field  Ea(r)  is  also  shown  for  comparison. 
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to  the  sheath  region  near  the  wall.  The  resulting  value  of  the  potential  at 
the  wall  relative  to  the  axis  is  V(r=R)a:!-6.6  V.  This  value  can  be  compared 
with  the  value,  Va,  given  by  Eq.  (22a)  in  Ref.  13  for  the  quasineutral  inertia 
limited  theory  of  the  positive  column.  For  the  discharge  parameters 
discussed  here,  Va=-9.5  V.  The  discrepancy  between  these  two  values  is  a 
result  of  the  assumptions  of  quasineutrality  and  inertia  limited  flow  used  in 
deriving  the  expression  for  Va. 

Figure  3  is  a  plot  of  the  radial  charged  particle  current  density 
Jr(r)=err(r).  The  maximum  in  this  curve  is  a  consequence  of  the  cylindrical 
geometry  and  the  fact  that  the  particle  production  rate  vine(r)  decreases 
with  increasing  r.  However,  the  total  radial  charged  particle  current  per 
unit  length  at  radius  r,  27crJr(r),  is  a  monotonically  Increasing  function  of  r, 
as  required. 

For  comparison,  the  results  of  calculations  for  R=0.5  cm  and  Po=1.0  Torr 
are  given  in  Figs.  4-6.  This  leads  to  a  value  cm-Torr,  which  is  10 

times  larger  than  the  value  of  Rpo  for  the  previous  case.  For  these 
conditions,  the  measurements  of  Lelby  and  Rogers  imply  Te=45,000  oK, 
leading  to  a  required  net  charged  particle  density  on  the  axis  of  An(r=!0)=4.9 
X1013  m'3.  The  resulting  ionization  frequency  is  Vi=6.67  X  lO^  s*l.  Figure  4 
shows  the  plots  of  nifr).  ne(r),  and  An(r)  for  this  situation.  As  in  the 
previous  case,  the  curves  for  ni(r)  and  n^Cr)  are  close  to  the  shape  of  the 
zero  order  Bessel  function  Jo(2.4r/R)  for  the  bulk  of  the  plasma,  with  a 
significant  deviation  occurring,  especially  for  ni(r),  only  near  the  tube  waU. 
This  deviation  is  again  indicative  of  the  development  of  a  sheath  region,  as  is 
also  shown  by  the  curve  for  An(r)  A  comparison  between  Fig.  1  and  Fig.  4 
shows  that  the  resulting  sheath  region  is  larger  for  smaller  values  of  Rpo.  In 
addition,  this  comparison  also  shows  that  the  values  for  An(r)  Increase  with 
decreasing  Rpo,  with  nearly  an  order  of  magnitude  difference  in  An(r) 
occurring  between  these  two  cases. 

Figure  5  presents  the  results  for  EJr(r)  and  V(r)  for  this  second  case, 
along  with  the  curve  for  Ea(r)  from  the  ambipolar  diffusion  theory.  Again, 

^he  curve  for  V(r)  demonstrates  the  smooth  transition  from  the  bulk  plasma 
’•egion  near  the  axis  to  the  sheath  region  near  the  wall.  For  this  situation, 
the  resulting  wall  potential  is  V(r=R)=-12.1  V,  as  compared  to  the 


13s.  A.  Self,  Phys.  Fluids  6.  1762  (1963). 


101 


n  (r)  (10^®  m"®) 


1.0 


1.30 


Fig.  4.  Numerical  calculations  of  the  ion  density  ni(r),  the  electron  density 
ne(r),  and  the  net  charged  particle  density  An(r)  for  Rpo=0.5  cm-Torr. 
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Fig.  5.  Numerical  calculations  of  the  radial  electric  field  strength  E^(r)  and 
the  resulting  electric  potential  V(r)  for  Rpo=0,5  cm-Torr.  The  ambipolair 
radial  electric  field  Ea(r)  is  also  shown  for  comparison. 
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quasineutral  inertia  limited  value  Va=-16.5  V.  The  magnitude  of  the  radial 
electric  field  shown  in  this  figure  is  significantly  larger  than  the  electric 
field  in  Fig.  2.  This  is  a  consequence  of  the  higher  electron  temperature  for 
the  situation  depicted  in  Fig.  5,  which  leads  to  enhanced  electron  diffusion 
to  the  wall.  This  results  in  the  generation  of  a  larger  space  charge  induced 
radial  electric  field  in  order  to  retard  the  electrons  and  accelerate  the  ions 
toward  the  wall.  As  a  result,  the  point  Vdi=Vri  is  closer  to  the  axis  in  Fig.  5 
than  in  Fig.  2,  since  the  ion  drift  velocity  at  a  given  relative  radial  distance 
from  the  axis  is  larger  for  a  larger  electric  field. 

Figure  6  shows  Jr(r)  for  the  situation  where  Rpo=0.5  cm-Torr.  The  value 
of  Jr(r)  for  this  case  is  more  than  an  order  of  magnitude  larger  than  those 
shown  in  Fig.  3. 

From  the  above  results,  it  is  possible  to  investigate  the  applicability  of  the 
Bohm  criterion,  Vdi^(kTe/mi)i/2^  for  these  discharge  conditions, For  both 
discharge  conditions  discussed  above,  Vdi(r)<(kTe/mi)i/2  for  all  values  of  r. 
Thus,  the  Bohm  criterion  is  never  met.  even  though  the  results  presented  in 
Figs.  1*6  clearly  show  the  development  of  a  sheath  region  near  the  discharge 
tube  wall.  The  formation  of  a  sheath  region  even  in  the  absence  of  the  Bohm 
criterion  being  fulfilled,  has  been  discussed  previously  by  Ingold.® 

Finally,  the  self-consistency  of  the  theory  in  terms  of  the  assumptions  and 
boundsuy  conditions  used  can  be  discussed  in  view  of  the  calculations 
presented  above.  Assumption  (4)  in  Sec.  II  states  that  vi  is  small  with 
respect  to  Vmi  and  Vme-  As  can  be  seen  by  comparing  the  calculated  values 
for  vi  with  the  v£dues  given  for  Vmi  and  Vmc.  this  assumption  is  indeed 
satisfied  for  the  cases  presented.  The  importance  of  the  boundary  condition 
chosen  for  the  electron  density  at  the  wall  in  determining  the  results  of  the 
calculations  can  be  determined  by  comparing  the  tube  radii  used  above  with 
the  tube  radii  predicted  by  other  commonly  employed  boundary  conditions 
for  the  same  discharge  parameters.  In  addition  to  the  botmdary  condition  of 
Eq.  (5)  which  was  used  in  these  calculations,  another  commonly  used 
boundary  condition  is  ne(r=Fd=0.  From  Figs.  1  emd  3  it  is  clear  that  this 
condition  is  closely  approximated  at  r=R  in  these  figures.  Indeed, 
calculations  show  that  the  tube  radius  is  Increased  by  less  than  0.2%  by  the 


Bohm,  In  The  Characteristics  of  Electrical  Discharges  tn  Magnetic  Fields,  edited  by 
Guthrte  and  R.  K.  Wakerling  (McGraw-Hill.  New  YorlC  1949).  Chap.  3. 
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use  of  this  boundaiy  condition.  Therefore,  the  assumption  of  the  boundary 
condition  of  Eq.  (5)  is  consistent  with  the  results. 


IV.  CONCLUSIONS 

A  model  heis  been  presented  for  the  radial  distribution  of  the  ion  and 
electron  densities  and  electric  held  strength  in  the  cylindrical  positive 
column  of  a  dc  discharge  for  a  plasma  consisting  of  a  singly  charged  positive 
ion  species,  electrons,  and  neutrals.  Under  appropriate  assumptions,  the 
resulting  set  of  equations  involved  consists  of  the  particle  conservation 
equation,  the  momentum  conservation  equations,  and  Poisson's  equation.  A 
method  of  solution  for  this  set  of  equations  using  a  power  series  expansion 
has  been  outlined,  which  avoids  the  instability  problem  encountered  by 
previous  authors.  This  allows  for  the  inclusion  of  all  the  volume  force  terms 
in  the  conservation  of  momentum  equations  for  the  ions  and  electrons. 

The  results  of  calculations  based  on  this  model  shew  the  development  of 
the  bulk  plasma  region  near  the  discharge  tube  axis  and  the  sheath  region 
near  the  tube  wall.  There  is  no  clear  distinction  between  these  two  regions, 
with  a  smooth  transition  occurring  from  the  bulk  plasma  region  to  the 
sheath  region,  as  can  be  expected.  While  the  results  near  the  tube  axis  are 
consistent  with  the  ambipolar  model,  they  deviate  considerably  from  this 
model  in  the  sheath  region  near  the  wall.  However,  even  at  the  axis,  there 
is  a  net  space  charge  density  which  increases  for  decreasing  values  of 
The  electric  potential  at  the  wall  relative  to  the  axis  is  on  the  order  of  the 
value  predicted  by  the  quasineutral  inertia  limited  theory,  with  the  waU 
potential  increasing  for  decreasing  values  of  Rpo.  The  results  show  that  the 
development  of  a  sheath  region  can  occur  even  for  situations  in  which  the 
Bohm  criteria  is  never  met. 

The  importance  of  the  boundary  condition  chosen  for  the  electron 
density  at  the  discharge  tube  wall  was  investigated.  It  was  found  that  the 
choice  of  either  Eq.  (5)  or  n€(r=IU=0  as  the  boundary  condition  has  little 
effect  on  the  solution  for  the  cases  studied. 
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